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On  p  11-23  the  top  line  should  read: 

In  regard  to  the  determination  of  S.,  figures  5,  7,  and  8  show  that 
2.2  was  slightly  better  than  2.0  or  2.5. 


On  p  II  -  23,  the  expression  for  K,  should  read: 
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The  values  of  the  vibrational  relaxation  time  given  in  Table  I  p  II  -  36 
were  not  reduced  in  the  same  fashion  as  Blackman's  data.  When  the  data 
was  reduced  by  the  identical  method  shorter  times  were  obtained  which 
agree  quite  well  with  his  data.  The  following  changes  should  be  made  in 
Table  I: 

For  run  No.  2  the  shock  velocity  should  be  2.96  mm//isec. 

The  new  relaxation  times  for  1  atmosphere  of  pressure  in  shock  fixed 
coordinates  in  ^lsec  are  respectively  for  run  Nos.  1  through  5:  2.6,  12.1, 
4.3,  2.7,  and  .2. 

The  inferred  efficiencies  of  Ne  and  A  in  exciting  vibrationally  remain 
unchanged.  All  the  N^  -  Ne  interferograms  showed  two  curvatures  in  the 
relaxation  zone  inferring  a  slow  vibrational  relaxation  prior  to  dissociation. 
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GENERAL  INTRODUCTION 


During  high  altitude  flight  at  hypersonic  speeds,  the  shock  wave 
which  forms  upstream  of  an  entry  vehicle  will  have  a  finite  thickness 
which  can  be  an  appreciable  portion  of  the  shock  detachment  distance, 
depending  on  the  geometry  of  the  vehicle,  the  flight  speed  and  the  free 
stream  density.  Further,  due  to  the  low  ambient  gas  density,  and  the 
consequent  lower  number  of  intermolecular  collisions,  there  is  insuf¬ 
ficient  time  for  the  gas  to  reach  a  new  state  of  thermochemical  equi¬ 
librium  during  the  residence  of  the  particles,  and  hence  the  gas  flow 
in  the  immediate  vicinity  of  the  vehicle  is  a  non-equilibrium  viscous  flow. 

In  Earth's  low  altitude  flight  regime  (below  150,000  ft.  ),  in  first 
approximation,  one  can  uncouple  the  fluid  dynamic  processes  from  the 
chemical  kinetics.  Here,  one  first  merely  solves  the  well-known 
algebraic  relationships  for  the  conservation  of  mass,  momentum  and 
energy  (Rankine-Hugoniot  equations)  and  then  proceeds  to  calculate 
the  rate  of  change  of  gas  composition,  temperature,  density  and  pressure 
by  simply  integrating  a  set  of  chemical  kinetic  equations. 

However,  in  the  regime  of  interest  such  a  procedure  is  not  meaning¬ 
ful  since  one  is  dealing  with  coupled  physicochemical  phenomena.  In 
particular,  the  available  kinetic  energy  of  the  molecules  is  redistributed 
among  many  different  internal  modes  of  energy  storage  which  absorb  energy 
at  different  rates,  but  nevertheless  in  parallel  processes  and  not  in  series. 
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The  theoretical  approach  to  this  problem  therefore  requires  the 
formulation  and  numerical  solution  of  the  coupled  conservation  laws, 
suitable  for  the  description  of  non-equilibrium  chemically  reacting  fluid 
flow  in  which  large  gradients  of  velocity,  temperature  and  pressure  are 
present.  This  type  of  theoretical  model  requires  the  inclusion  of 
molecular  transport  processes,  relaxation  phenomena  and  chemical 
kinetics.  In  order  to  accomplish  this  task,  one  must  have  a  knowledge 
of  the  various  physicochemical  quantities  which  appear  as  variable 
coefficients  in  the  governing  differential  equations. 

The  experimental  approach  consists  of  the  development  of  a  shock 
tube  facility  in  which  accurate  relaxation  rate  data  and  chemical  kinetic 
data  may  be  obtained.  These  data  can  then  be  utilized  to  supplement  and 
verify  the  macroscopic  predictions  of  the  theoretical  investigation. 

The  following  two  sections  contain  the  results  of  Phase  I  of  a 
theoretical  and  experimental  study  of  shock  wave  structure  in  a  relaxing, 
dissociating  diatomic  gas.  The  theoretical  study  was  carried  out  by 
Scala  and  Talbot,  and  appears  as  Part  1.  The  experimental  investigation 
was  carried  out  by  Cary  and  appears  as  Part  II. 

It  should  be  notea  that  the  theory  is  quite  general  and  applies  to  an 
arbitrary  diatomic  gas,  while  the  experiments  were  performed  utilizing 
molecular  nitrogen,  one  of  the  primary  constituents  of  Earth's  atmosphere 
and  the  Martian  and  Venusian  atmospheres. 
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ABSTRACT 


In  an  earlier  paper  (Ref.  1)  the  authors  analyzed  the  structure  of  a 
shock  wave  in  a  gas  having  a  long  relaxation  time  for  its  internal  molecular 
(rotational)  degrees  of  freedom.  This  analysis  has  now  been  extended  to 
the  case  of  a  diatomic  gas  having  two  relaxation  effects,  one  for  the 
rotational  and  one  for  the  vibrational  degrees  of  freedom.  Departures 
from  thermodynamic  equilibrium  within  the  shock  wave  were  studied  for 
two  values  of  the  Mach  number,  M  =  4.  6  and  M  =  10.  The  number  of 

inter  molecular  collisions  required  for  adjustment  of  the  rotational  degrees 
of  freedom  was  kept  fixed,  but  the  number  of  collisions  required  for 
vibrational  adjustment  was  varied  parametrically  from  300  to  2000. 

The  relaxation  model  chosen  to  supplement  the  Navier-Stokes 
equations  was  a  generalization  of  that  used  in  reference  1,  namely,  both 
the  rotational  and  vibrational  temperatures  Tj.  and  Ty  change  through  the 
shock  wave  at  a  rate  which  is  proportional  to  the  difference  between  them¬ 
selves  and  the  other  two  temperatures,  and  inversely  proportional  to  the 
number  of  intermolecular  collisions.  In  this  model,  there  is  direct 
coupling  between  Tr  and  Ty,  as  well  as  indirect  coupling  through  the 
translational  temperature  Tj. 

The  calculations  provide  a  more  realistic  estimate  for  the  distribution 
of  translational  temperature  T^  through  a  normal  shock  wave  than  is  obtained 
with  the  neglect  of  relaxation  phenomena.  Accurate  estimates  of  T^  are  of 
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importance  in  determining  the  onset  and  initial  rates  of  chemical  reactions 
in  strong  shocks. 

The  calculations  also  give  a  more  realistic  value  of  the  effective 
thickness  of  the  shock  transition.  This  information  is  of  importance  in 
the  problem  of  the  merging  of  the  shock  wave  with  the  shock  layer  at  the 
nose  of  a  blunt  body  in  high  altitude  hypersonic  flow. 
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INTRODUCTION 


The  shock  wave  structure  problem  has  been  of  interest  to  aero- 
dynamicists  for  many  years,  and  this  interest  has  been  heightened 
recently  with  the  advent  of  hypersonic  flight.  A  considerable  body  of 
theoretical  and  experimental  work  therefore  exists  which  has  been 
reviewed  recently  by  Sherman  and  Talbot  (Ref.  2),  Talbot  (Ref.  3)  and 
Losev  and  Osipov  (Ref.  4).  The  bulk  of  the  work  is  confined,  however, 
to  the  study  of  shock  wave  structure  in  non-reacting  monatomic  gases 
which  therefore  gives  little  insight  into  the  problem  of  the  shock  structure 
in  a  real  gas  at  hypersonic  speeds. 

In  a  diatomic  gas,  one  must  consider  the  contributions  of  the 
internal  degrees  of  freedom  (e.  g.  rotational,  vibrational)  as  well  as 
translational  in  studying  the  approach  to  thermal  equilibrium  through 
a  shock  wave.  Further,  the  effects  of  finite  chemical  kinetics  (dissociation 
and  ionization)  which  begin  to  appear  at  the  Mach  numbers  associated  with 
hypersonic  flight  are  coupled  with  the  other  physicochemical  processes, 
and  thus  the  problem  requires  the  simultaneous  solution  of  the  coupled 
fluid  dynamic  and  chemical  kinetic  equations.  This  is  particularly  true 
at  the  highest  flight  speeds  (highest  non-equilibrium  temperatures)  where 
the  relaxation  times  for  the  internal  processes  and  the  chemical  processes 
are  comparable. 

The  earliest  theoretical  studies  of  shock  waves  were  carried  out 
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for  a  non-reacting  perfect  fluid,  utilizing  the  Navier -Stokes  equations.  For 
example,  Rankine  (Ref.  5)  solved  the  problem  for  the  case  of  zero  viscosity 
and  constant  thermal  conductivity.  The  classical  Navier-Stokes  treatment 
is  due  to  Becker  (Ref.  6)  who  obtained  a  solution  when  the  viscosity  and 
thermal  conductivity  were  assumed  independent  of  the  gas  temperature. 
Liater,  Thomas  (Ref.  7)  extended  the  work  of  Becker  and  showed  that  if 
the  viscosity  of  the  gas  is  proportional  to  the  square  root  of  the  absolute 
temperature,  then  a  strong  shock  wave  solution  results  in  a  non-vanishing 
shock  thickness. 

Recently,  Grad  (Ref.  8)  and  Gilbarg  and  Paolucci  (Ref.  9)  obtained 
solutions  for  the  Navier-Stokes  equations  utilizing  more  realistic  values 
of  the  specific  heats  ratio  and  the  Prandtl  number,  where  the  viscosity 
was  assumed  to  follow  a  power  law  dependence  on  the  absolute  temperature. 
Their  results,  svhon  compared  with  available  experimental  data,  showed 
the  importance  of  selecting  correct  transport  coefficients  in  calculating 
shock  wave  thickness  and  have  given  greater  credence  to  the  Navier- 
Stokes  model  for  shock  wave  structure,  at  least  for  weak  shocks. 

With  regard  to  internal  relaxation  phenomena  associated  with  the 
internal  degrees  of  freedom  of  the  molecules,  there  is  the  early  work 
o*"  Bcthc  and  Teller  (Ref.  10),  Broer  (Ref.  11)  and  Bighthill  (Ref.  12)  in 
which,  however,  the  shock  transition  is  considered  to  occur  in  two 
distinct  steps. 
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More  recently,  studies  have  been  published  on  relaxing  shock  waves 
in  which  the  relaxation  phenomena  is  taken  into  consideration  throughout 
the  entire  shock  transition.  The  work  of  Broer  and  van  der  Bergen  (Ref. 
13)  is  valid  only  for  very  weak  shocks,  however,  the  treatment  of  Talbot 
and  Scala  (Ref.  1)  based  on  the  Navier-Stokes  equations  and  an  auxilliary 
relaxation  equation  seems  to  be  phenomenologically  correct  for  strong 
shock  waves. 

In  reference  1,  the  authors  analyzed  the  structure  of  a  shock  wave 
in  a  diatomic  gas  in  which  the  rotational  degrees  of  freedom  were  assumed 
to  adjust  more  slowly  than  the  translational  degrees  of  freedom,  and  the 
vibrational  degrees  of  freedom  were  unexcited.  The  problem  was  treated 
according  to  both  the  bulk  viscosity  and  the  relaxation  equation  points  of 
view.  It  was  shown  that  the  bulk  viscosity  model  was  limited  to  such 
small  values  of  the  relaxation  time  as  to  preclude  its  application  to  most 
cases  of  interest,  whereas  the  relaxation  equation  model  yields  results 
which  seem  reasonable  for  the  full  range  of  relaxation  times  and  Mach 
numbers  treated,  0  <  N  <  100  and  I  <  M  <  10. 

In  the  present  paper,  the  analysis  has  been  extended  to  the  case 
of  a  diatomic  gas  having  two  relaxation  effects,  one  for  the  rotational 
and  one  for  the  vibrational  degrees  of  freedom.  Departures  from  thermo¬ 
dynamic  equilibrium  within  tlie  shock  wave  were  studied  for  two  values 
of  Mach  number,  M  =  4.  6  and  10.  The  number  of  intermolecular 
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collisions  required  for  adjustment  of  the  rotational  degrees  of  freedom 
was  kept  fixed,  but  the  number  of  collisions  required  for  vibrational 
equilibration  was  varied  parametrically  from  300  to  2000  in  order  to 
study  the  coupling  between  the  translational,  rotational  and  vibrational 
temperatures. 

For  the  sake  of  generality,  the  fundamental  equations 
have  been  derived  for  the  one-dimensional  flow  of  a  dissociating 
diatomic  gas  having  rotational  and  vibrational  relaxation.  In  the 
present  work,  however,  solutions  including  finite  dissociation  rates 
were  not  investigated. 
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NOMENCLATURE 


Si  f  Id  f  c  f  d  t  6 


C 


e 


constants 

dimensionless  specific  heat  at  constant 
volume 

mean  molecular  speed 

mass  fraction  of  atomic  species 

specific  heat  at  constant  volume  | 

binary  diffusion  coefficient 

thermal  diffusion  coefficient  of  species  i 

internal  energy 
energy  of  formation 
enthalpy,  including  chemical 

heat  of  formation  | 

diffusion  flux 

I 

thermal  diffusion  ratio 

mean  free  path 
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i 


I 


m 

M. 

1 

M 

N 

P 

P 

Q 


R 


Ri 


u 


V 


w 


n 


reference  length  based  on  viscosity 


integration  constant  in  continuity  equation 


molecular  weight  of  species  i 


mean  molecular  weight 


number  of  intermolecula r  collisions 


static  pressure 


integration  constant  in  momentum  equation 


integration  constant  in  energy  equation 


dimensionless  density 


universal  gas  constant 


gas  constant  of  species  i 


temperature 


velocity 


diffusion  velocity 


chemical  source  term 


total  molar  density 
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X 


physical  spatial  variable 


y 

Z 

a ,  fi,  y,  S,  ( 

K 

X 

A 

P 

6 

T 

4> 

a 

Subsc  ripts 


dimensionless  spatial  variable 
integration  constant  in  diffusion  equation 


constants 


bulk  viscosity  coefficient 
thermal  conductivity 


viscosity 


density 

dimensionless  group 

dimensionless  translational  temperature 
dimensionless  rotational  temperature 
dimensionless  vibrational  temperature 


dimensionless  pressure 


A 


atoms 


species  1 


molecule  s 


rotational 


translational 


vibrational 


upstream 


mode  of  energy  storage 


downstream 


reference  state 
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BASIC  EQUATIONS 


It  is  assumed  that  the  total  thermal  energy  of  the  gas  is  composed 
of  translational,  rotational  and  vibrational  modes,  and  hence  upon  char¬ 
acterizing  these  different  modes  by  temperatures  Tj,  T^.  and  and  specific 
heats  and  ,  the  internal  energy  of  the  gas  is 

e  =  ?  C.  e.  ;  de.  =  2  C  dT.  { 1) 

111  1  n  V.  1 

1  n 


in  which  case  the  ordinary  specific  heat  at  constant  volume  is: 

C  =  C  +  C  +  C 

V.  V.  V.  V.  (2) 

1111  '  ' 
t  r  V 

For  convenience  it  is  assumed  that  all  specific  heats  are  constant. 

For  steady,  one-dimensional  flow,  the  following  conservation  laws 
apply.  The  conservation  of  mass  is  given  by: 


(  p  u)  =  0 
dx 

which  upon  integration  becomes: 

p  u  =  m  (a  constant) 

Tlie  conservation  of  momentum  is  given  by: 

p  u  du  =  -  ^  +d  r{4p+K)du 
d.N  dx  dx  L  3  dx 

whicli  when  integrated  becomes: 

mu  +  p-{4^+«)du  =  P  (a  constant) 

3  ^ 

The  conservation  of  species  i  is  given  by: 


(3) 

(4) 

(5) 

(6) 


_d  (  puC.  +  )  =  w.  (7) 

dx 


1 


which  when  integrated  becomes: 


mq  +  jj  -  $  i  =  Z  (a  constant) 

where: 

d  $  .  =  w.dx 

and  the  one -dimensional  diffusion  flux  vector  for  a  binary  mixture  may  be 


written  (Ref.  14): 
j  =  -  P  ^ 


. .  j  dCi  +  n  M  -  Mi  )  C  .  (]  1 

L  Mi 

■Cq-  d  In  I 


-f  K'p 


where  the  thermal  diffusion  ratio  is  given  by  (Ref.  15)  : 


(10) 


KT=__jf_  dJ  (") 

n^MiMj  "^ij 

The  energy  equation  may  be  wrHtcn: 


p  u  dh  = 
dx 

_d_  ( 

dx 

A,  dT,  + 
dx 

^r 

dT  + 

■1  _r 

dx 

dT,,  ) 
dx 

-JL(Z 

d.x  i 

ji^i 

)  +  (  4 
3 

+  K  ) 

(  ^)^ 
dx 

(12) 

+  u  + 

d 

■  RT^ 

z 

Z 

dx 

dx 

n 

i 

j'i 

M.-©..  J 

1  ij 

Upon  introducing  the  momentum  equation  (5)  and  integrating,  this  becomes 

2mh  +  mu  -  2  A  ^  dT^  -  2A^d1^  -  2A^dl^ 

Ox  dx  dx 


2  Xj^hi  -  2u  (  4  ^  +  K  )  ^ 


-  2 


1 

RT, 

n 


n.D; 


i 


dx 


(V.  -  V.) 


(13) 


Q  (a  constant) 
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The  thermal  conductivities  A  ^ ,  A  j.  and  A  ^  corresponding  to  the 
transport  of  translational,  rotational  and  vibrational  modes  of  energy 
are  evaluated  according  to  the  method  of  Eucken,  (See  Ref.  15). 

For  a  binary  mixture  of  atoms  and  molecules,  the  Dufour  term 
may  be  expanded  to  yield; 


-  2  RT 


n  i  j  /“i  M.  A  . 


1  ij 


Noting  that  the  enthalpy  may  be  written; 


(14) 


h  =  e  f  p 

P 


(15) 


and  that  the  internal  energy  of  the  i  ili  species  is  given  by: 
T 

«i  =  fc  dTi  +  Aef° 

V  1  1 

o 

The  energy  equation  becomes; 

2m  (C  T  +  C  T  +  C„  T.  )  +  2pu 

V.,  t  V.,  r  V  V  ' 

Mj.  l^v 

+  mu^  -2A,  dT,  2X  dT  -2  A  dT 

t  t-  r  r  V  V 


dx 


dx 


dx 


2u(4p  +  K  )  d^  +  2mCAhf° 


+  2mC 


.  [<c„ 

r  (c  +  R 

i 


dx 
-  C 


)  T  -  C  T  -  C 
V  X  €  t  V,  ,  r  \ 

Mj. 


(16) 


(IV) 
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Vx  f  V 
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.  o  Pl^T 
T..  -(  Ahr  +  — t:: T-,  _  . 


] 


=  Q 
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The  thermodynamic  pressure  is  identified  with  the  translational 
temperature  and  hence  the  equation  of  state  may  be  written: 
p  =  P  RT^ 
or  equivalently; 

p  “  "  V  +  =a' ’'t  (w 

The  relaxation  equations,  following  Wang  Chang  and  Uhlenbeck, 
(Ref.  16).  but  extending  the  concept  to  a  three-temperature  model,  may 
be  put  into  a  symmetric  form  (Refs.  17  and  18)  which  may  be  written: 


udT  =T-T+T-T  (21) 

V  t  V  r  V 


where  the  relaxation  time  (8  may  be  related  to  the  mean  free  path  I* 

the  mean  molecular  speed  C  and  to  the  number  of  intermolecular  collisions 

for  relaxation  N  through 

/3  =  =  J[_  N  P  (22) 

C  4  p 

Also,  in  the  derivation  of  the  equations,  it  has  been  assumed  throughout 
that  the  temperature  of  the  atoms  is  equal  to  the  translational  temperature 
of  the  molecules. 
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SOLUTION  OF  GOVERNING  EQUATIONS 


Equations  A.  7  through  A.  12  given  in  the  Appendix  were 
integrated  numerically  in  phase  space  by  means  of  an  electronic  analog 
computer  for  two  values  of  the  Mach  number,  4.  6  and  10.  In  order 
to  examine  the  coupling  between  the  three  temperatures,  t  ,  and  0, 
the  number  of  collisions  required  for  vibrational  excitation  was  varied 
in  several  steps  between  300  and  2000  intermolecular  collisions.  The 
results  of  the  phase  space  integrations  are  suown  in  Figures  1  through  6. 

It  is  seen  that  with  the  form  of  relaxation  equations  given 
in  equations  A.  12,  the  rotational  temperature  starts  to  rise  almost  imme¬ 
diately  as  the  translational  temperature  rises,  however,  at  a  slower  rate 
of  excitation.  The  rotational  temperature  very  closely  approaches 
the  translational  temperature  r  but  never  exceeds  it  which  is  a  satisfactory 
phenomenological  result.  The  vibrational  temperature  ^  lags  furthest 
behind  and  its  effect  is  felt  at  about  the  same  time  both  the  rotational  and 
translational  temperature  go  through  a  maximum.  For  the  smallest 

number  of  collisions  considered  (N  =  300),i/f  actually  starts  rising  before 

vr 

the  rotational  temperature  reaches  a  maximum. 

In  general,  then,  one  observes  that  r  and  <f>  are  strongly 
coupled  while  the  coupling  between  ^  and  both  r  and  <t>  is  weaker. 

The  slower  rise  in  the  vibrational  temperature  acts  to  produce  a  relaxation 
"tail"  in  both  the  translational  and  rotational  temperatures  which  overshoot 
their  final  equilibrium  values.  On  the  other  hand,  if/  behaves  in  a 


1 


monotonic  fashion  for  all  values 
collision  s. 


of  N 

vr 


in  a  range  between  300  and  2000 


Although  an  increase  in  Mach  number  acts  to  change  the 
downstream  boundary  conditions  relative  to  the  upstream  boundary 
conditions,  qualitatively,  it  was  found  that  there  is  no  significant 
change  in  the  structure  of  a  normal  shock  in  the  range  of  Mach  numbers 
between  4.  6  and  10.  The  numerical  solutions  shown  in  Figures  1 
through  6  were  transformed  back  into  the  spatial  plane  by  quadrature 
by  means  of  an  IBM  7090  digital  computer.  In  the  spatial  plane,  the 
center  of  the  coordinate  system,  y  =  0,  was  arbitrarily  selected  to 
coincide  with  that  point  within  the  shock  wave  at  which  the  velocity 
takes  on  its  average  value,  i.  e.  ,  v  =  1/2  (v.^  +  v.,oe)-  Here  it  is  seen 
that  althovigh  the  thermodynamic  pressure  is  associated  with  the 
translational  temperature,  neither  the  pressure  a  nor  the  density  r 
exhibit  the  o\'ershoot  that  the  translational  temperature  e.xhibits.  That 
is,  with  the  e.xception  of  r  and  0  all  physical  variables  e.xhibit  a  monotonic 
behavior  for  the  full  range  of  conditions  considered  here  (see  Figures  7 
through  12).  In  obtaining  these  solutions  it  has  been  assumed  that 
dissociation  does  not  occur. 
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SHOCK  WAVE  THICKNESS 

One  nomin5.l  measure  of  shock  thickness  is  based  on  the  maximum 
slops  of  one  of  the  physical  variables.  That  is,  for  a  given  property  f 
(which  may  be  taken  conveniently  as  T,  <f> ,  v,  O  or  r  ),  the  maximum 
slope  thickness  may  be  expressed  in  nondimensional  form  as 

K  _  I  ^max  ”  ^minl 

■  pdf/dV)max] 

The  maximum  and  minimum  values  of  the  vibrational  temperature,  gas 
velocity,  pressure  and  density,  (i.  c.  ,  v,  0  and  r)  occur  at  y  =  ^  ; 

however,  the  maximum  values  of  the  Iranslatioiial  arid  rotational  tempera¬ 
tures  (i.e.  T  and  <f))  will  usually  occur  at  some  intermediate  location, 

(unless  all  of  the  relaxation  times  arc  assumed  to  be  identically  /-ero). 

The  term  |  (df/dy),„jt,<  |  represents  the  absolute  value  of  the  maximum 
gradient  of  f,  regardless  of  where  within  the  shock  wave  it  occurs. 

Figs.  13  and  I -I  show  that  the  maximum  slop  thicknesses  based  on  t. 
(f),  V  ,  <r  and  r  are  essentially  independent  of  Ny  .for  a  given  Mach  number 
and  Nr,  whereas  (Ay)4>  increases  directly  and  almost  linearly  with  N^  ,  • 
Since  Nr  was  hold  fixed  while  Ny^  was  varied  parametrically,  it  is  not 
surprising  lhat(Ay)4>  hchaves  differently  for  the  three-temperature  model 
treated  here  than  in  the  two  temperature  model  treated  earlier  in  reference 
1.  That  is,  it  is  found  here  that  in  spite  o'"  the  direct  coupling  between  the 
rotational  and  vibrational  tenr peraturcs  in  the  two  relaxation  equations,  0 
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primarily  follows  r  and  goes  through  a  maximum.  Thus,  the  primary 
effect  of  an  increase  in  is  to  increase  the  maximum-slope  thickness 

based  on  <A.  and  to  simultaneously  add  a  long  relaxation  "tail"  to  r  and 
<ft,  which  cannot  appreciably  affect  the  value  of(Ay)^  . 

In  order  to  obtain  an  appreciation  of  the  implications  of  these  new 
results  upon  the  thickness  of  a  shock  wave  in  Earth's  atmosphere,  the 
values  of  the  dimensionless  shock  thickness  (Ay)^  shown  in  Figure  13  were 
transformed  into  the  physical  thicknesses  (A^^)^  shown  in  Table  II.  Simulated 
initial  conditions  corresponding  to  a  range  of  altitudes  between  ZOO,  000  and 
400,000  ft.  were  utiliiied,  and  it  was  found  that  at  a  Mach  number  of  10, 
and  =  2000  collisions,  the  thickness  of  the  shock  wave  corresponding 

to  equilibration  of  the  translational,  rotational  and  vibrational  degrees  of 
freedom  varied  by  five  orders  of  magnitude,  i.  e.  10  ^  <  (Ax)^  <  10  ft. 
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CONCLUSIONS 


A  new  theoretical  model  for  shock  wave  structure  in  a  diatomic  gas 
with  rotational  and  vibrational  excitation  was  developed  here,  and  numerical 
integrations  have  been  carried  oat  employing  high  speed  electronic  computing 
equipment.  The  model  which  is  based  on  the  use  of  symmetric  relaxation 
equations  yields  results  which  seem  reasonable  for  the  full  range  of 
relaxation  times  and  Mach  numbers  treated. 

More  realistic  estimates  are  provided  for  the  distribution  of  trans¬ 
lational  temperature  through  a  normal  shock  wave  than  would  have  been 
obtained  with  the  neglect  of  relaxation  phenomena.  These  results  will  be 
of  use  in  further  theoretical  studies  of  the  rates  of  chemical  reactions  (e.  g. 
thermal  dissociation)  in  strong  shocks. 

The  new  calculations  give  more  realistic  values  for  the  effective 
thickness  of  the  shock  transition,  which  will  yield  better  criteria  for  the 
determination  of  the  condition  of  "viscous  merging"  at  the  nose  of  a  blunt 
body  in  high  altitude  hypersonic  flow.  Examination  of  the  concept  of 
nominal  shock  wave  thickness  according  to  the  symmetric  relaxation 
model  shows  that  values  based  on  the  maximum  gradients  of  velocity, 
pressure,  density,  translational  temperature  and  rotational  temperature 
are  relatively  insensitive  to  the  magnitude  of  the  vibrational  relaxation  time, 
whereas  the  thickness  which  is  based  on  the  maximum  gradient  of  the 
vibrational  temperature  is  directly  proportional  to  the  number  of  intermolecular 
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collisions  required  for  vibrational  equilibration. 

It  can  also  be  inferred  from  these  calculations  that  in  the  initial 
steep  part  of  the  shock  wave,  direct  coupling  between  Tj.  and  is  not 
of  importance,  and  that  in  general  these  coupling  effects  affect  the 
temperature  distribution  only  in  the  relax.ation  "tail",  at  least  in  the 
range  of  parameters  and  flow  conditions  investigated.  The  coupling 
between  and  is  not  as  strong  as  between  and  Tj.,  but  does  exist 
and  depends  critically  on  the  magnitude  of  . 
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APPENDIX 


We  will  now  derive  the  non-dimensional  form  of  the  governing 
equations.  Following  Grad  (Ref.  8  ),  we  introduce  the  following  dimen¬ 
sionless  variables 


V 


mu 

P 


PP 

r  = - 

2 

m 


r 


<(>  = 


m^R 


P 


_M 

2 


mQ 


^  = 


m''R.,T 


M  *  V 


(A.  1) 


Let  us  introduce  the  transport  lengths 


m 


m 


P 


^t 

mRM 


(A.  2) 


1-19 


^  mRj^ 

L  = _ 

Now,  let  us  put  the  transport  property  equations  in  the  form: 


•*AM  = 


=  d(Tj' 


(A.  3) 


where  a,  b,  c,  d,  e,  a  ,  y ,  8  and  £  are  assumed  to  be  constants 

and  the  a  in  the  exponent  of  the  viscosity  law  is  not  to  be  confused  with 
the  a  appearing  in  equation  A.  1 , 

Upon  selecting  a  reference  temperature  we  can  introduce 

the  non-dimensional  spatial  variable: 
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y 


(A.  4) 


X 


where 


(A.  5) 


so  that  in  terms  of  the  new  variables  the  governing  equations  become  for 


the  case  of  difficult  energy  exchange  with  k  =  0 
mass: 

rv  =  1  (A,  6) 

momentum: 

4 
3 


a  dv 

f  V  — -- 

dy 


+  f(l  +  C.  ) 


(A.  7) 


species: 


Ja  =  -  ""^^A  -  ‘^A-oo^ 


(A.  8) 
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d/ 
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In  the  abov'e  eq'.iations, 
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The  dimensionless  ratios  6^,  6^  and  6^  are  defined  by: 


6. 


6  - 


,  sS  * 

o 


(A.  14) 


^  '  Lu"'  {t'‘)^-« 


It  may  readily  be  verified  that  the  ratios  6^,  Bj.  and  B^  have  the 
behavior  of  reciprocal  Prandtl  numbers,  whereas,  B  has  the  behavior 

a 

of  the  reciprocal  of  the  Schmidt  number. 

The  boundary  conditions  corresponding  to  the  Rankine-Hugoniot 
jump  relations  in  non-dimensional  form  may  be  obtained  by  setting  all 
the  derivatives  with  respect  to  y  equal  to  zero  in  equations  (A,  7)  through 
(A.  12). 
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TABLE  I.  CONSTANTS 


oc  (exponent)  =  V  =  5  =  <  =  1/2 

/3  (exponent)  =  3/2 
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rr 
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N  , 

rth 
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N 

vT 

= 

300, 500, 2000 

N  A 

VA 

= 

4000 

a  (M  =  4.6) 

= 

1.  238 

Q  (M  =  10) 

1.054 
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sed  on  maximum  slop  thickness  of  the  vibrational  temperature. 
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SHOCK  WAVE  PROFILES  IN  PHASE  SPACE  WITH 
ROTATIONAL  AND  VIBRATIONAL  RELAXATION 

0.18  r 


FIG.  6 


SHOCK  WAVE  STRUCTURE  IN  THE  SPATIAL  PLANE  WITH 
ROTATIONAL  AND  VIBRATIONAL  RELAXATION 
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INTRODUCTION 


This  investigation  was  undertaken  with  two  primary  objectives.  The  first 
objective  was  to  see  if  chemical  reaction  rates  can  be  determined  to  about 
the  same  accuracy  behind  combustion  driven  shock  waves  as  that  previously 
attained^’  ^  behind  shock  waves  generated  by  high  pressure  hydrogen  or 
helium.  The  second  was  to  obtain  chemical  rate  data  for  nitrogen.  The  following 
reaction  has  been  studied: 

+  M  ~  N  +  N  +  M 
where  M  =  N^,  N,  Ne,  or  A. 

The  above  reaction  becomes  of  increasing  importance  in  the  air  system 
for  re-entry  Mach  numbers  of  25  and  upwards  at  an  altitude  of  100,  000  feet 
and  above. 

The  earliest  experimental  study  of  nitrogen  dissociation  at  high  tempera¬ 
tures  was  carried  out  by  Eckerman'^.  A  hypervelocity  gun  was  used  to  fire 
small  spheres  into  N^  -  Xe  mixtures  at  high  velocities.  Shadowgraphs  were 
recorded  by  a  rotating  mirror  camera.  These  photographs  allowed  a  determina¬ 
tion  of  the  distance  from  the  front  of  the  sphere  to  its  bow  shock  wave. 

Eckerman  has  developed  a  method  for  relating  the  shock  detachment  distance 
to  the  bimolecular  dissociation  rates  of  a  diatomic  gas. 

In  1959,  Byron^  reported  interferometric  measurements  of  the  time  to 

reach  50%  of  equilibrium  dissociation  behind  shockwaves  in  nitrogen-argon 

mixtures.  His  technique  was  identical  to  that  used  previously  for  a  study  of 

.  .  2 

oxygen  dissociation. 


Allen  and  Camm  have  made  a  measurement  of  the  recombination  rate  of 
nitrogen  with  the  nitrogen  atom  and  molecule  as  third  bodies.  The  recombina¬ 
tion  rate  can  be  related  to  the  rate  of  decay  of  the  intensity  of  non-equilibrium 
radiation  behind  the  shock  front.  Hammerling"^  has  developed  a  theoretical 

model  for  the  origin  of  this  non-equilibrium  radiation. 

8  9  10 

At  low  temperatures,  Herron  ,  Wentinck  ,  and  Harteck  have  obtained 
recombination  rates.  Herron  measured  the  decrease  of  nitrogen  atom  concen¬ 
tration  as  a  function  of  distance  from  a  microwave  discharge  over  a  temperature 
range  of  195°K-453°K.  The  atom  concentrations  were  measured  by  "titration" 
with  NO^  and  determination  of  the  resulting  NO  by  mass  spectrometry.  Wentinck 
measured  the  recombination  rate  at  room  temperature  in  a  non-flow  system 
by  using  an  e  lectrodeless  discharge  to  partially  dissociate  the  nitrogen  in  a 
spherical  pyrex  flask  at  room  temperature.  A  platinum  resistance  thermometer 
was  centered  in  the  flask.  Rate  of  change  of  the  wire  resistance  could  be  cor¬ 
related  with  the  rate  of  disappearance  of  nitrogen  atoms.  Harteck  carried  out 
a  measurement  of  the  atom  recombination  rate  at  room  temperature  employing 
a  discharge  and  following  the  radiative  decay  of  the  yellow  nitrogen  afterglow. 

To  determine  the  recombination  rate,  it  is  necessary  to  measure  the  absolute 
atom  concentration  at  one  point.  Harteck  accomplished  this  by  "titration" 
with  NO. 

The  results  obtained  in  this  experiment  will  be  compared  with  those  just 
obtained  in  the  final  section. 
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THEORETICAL  CONSIDERATIONS 


If  thermal  and  mass  diffusion  are  not  large  throughout  the  dissociation 
zone,  and  moreover,  if  the  vibrational  degrees  equilibrate  before  appreciable 
dissociation  occurs,  then  it  is  possible  to  measure  inferred  specific  rate 
constants. 

The  equation  adopted  for  the  net  rate  of  decrease  in  the  concentration  of 
nitrogen  molecules  in  a  mixture  of  a  rare  gas  and  nitrogen  is 

Hi— 4  =  -  ■'d,  M 

2  3  2 

+  k  [n  ]  [n]  +k  [n]  +k  [n]  [Rare]  (1) 

1  *’2  *■3 

The  number  of  unknown  rate  constants  in  (1)  can  be  reduced  from  six  to  three, 
if  the  usual  assumption  is  invoked  that  at  equilibrium 

kd.  kd  kd, 

-r -  =  -r -  =  -r -  =  K 

kri  kr^  kr3 

where  K  is  the  equilibrium  constant.  Eq.  (1)  must  be  expressed  in  terms  of 
the  degree  of  dissociation  and  the  density  in  the  reaction  zone  behind  the 
shock  front  for  comparison  with  interferometric  measurements.  Mathews^ 
has  given  a  form  of  Eq.  (1)  which,  with  the  conservation  equations  of  mass, 
energy,  and  momentum,  the  shock  velocity,  and  the  equation  of  state,  yields 
the  variation  of  the  density  of  a  dissociating  diatomic  gas  as  a  function  of 
distance  behind  the  shock  front.  This  formulation  can  be  readily  extended  to 
include  mixtures  of  a  rare  and  a  diatomic  gas.  The  Appendix  gives  the  equations 


1 
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that  were  used  in  this  study  as  well  as  the  IBM  1620  program.  Since  the 
specific  refractivity  of  the  nitrogen  atom  differs  from  that  of  the  molecule 
the  program  was  set  up  to  generate  the  index  of  refraction  as  a  function  of 
distance  behind  the  shock  front. 

The  equilibrium  constant  has  to  be  evaluated  over  the  temperature  range 
of  these  experiments  (6,  000°K  -  8,  500°K).  Treanor's^^  values  are  tabulated 
in  1,000°K  steps  up  to  8,  000°K.  Baulknight's  tables  range  up  to  15,  000°K 

in  smaller  steps;  hence,  they  were  used  throughout. 
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Evans  has  shown  that  the  density  in  the  reaction  zone  can  be  expressed 
as  a  function  of  the  shock  velocity,  initial  gas  state  ahead  of  the  shock,  and 
a  dimensionless  parameter  /9  defined  by 


the  enthalpy  of  an  ideal  gas  at  0  K 


H 

R 


the  enthalpy  of  the  gas  at  temperature  T 

c  &  1  c  &  1 

the  universal  gas  constant  in  - —  if  H  is  expressed  in 


mol°K 


mol 
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Fickett  has  tabulated  the  values  oi/Slor  the  nitrogen  molecule  up  to  12,000°K 
which  is  adequate  for  the  temperature  range  carried  in  these  experiments. 

P  for  atomic  nitrogen  was  calculated  from  5,  000°K  to  12,  000°K  using 
the  following  relation: 

“  F.  ' 

I  i. 


n 


P  =  5/2  +  ^ 

n  ' 


n  =  1 


’n  k  T 
o 


e  k  T 
o 


4 

I 

n  =  1 


k  T 
o 


(2) 
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statistical  weight  of  the  nth  energy  level  of  the  nitrogen  atom 
the  energy  of  the  nth  energy  level  of  the  nitrogen  atom 
the  Boltzmann  factor 


T  =  the  absolute  temperature 

Collision  rate  expressions  were  used  for  i  and  k^  given  in 

Eq.  (1).  These  are: 


^^^«1  /  -R  N*'*"  .1/z/  D  ..„r  D  1 

(s^  +  i)  \mw^  j  ^  (rt]  p  [■  rtJ 


r  {S2  +  1) 


6  »r  R 
M  W,. 


N^-N 


ZffR 

MW_  •  MW 


N  -  Rare 
Z 


Rare 


3  ■r(S3+l) 


where  the  rate  constants  are  expressed  in 

^Z’  ^3  ~  collision  efficiency  factors 

N  =  Avogadro's  number 
o 

=  diameter  of  the  molecule 
5^  =  diameter  of  the  nitrogen  atom  taken  equal  to 
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the  average  of  the  squares  of  the  diameters  of  the  molecule  and 


rare  gas  atom 

D  =  dissociation  energy 

R  =  universal  gas  constant 

T  =  absolute  temperature 

Sj ,  =  one-half  of  the  square  terms  available  in  internal  energy 

of  the  molecule,  nitrogen  atom,  and  rare  gas  atom  respectively. 

The  above  expressions  are  based  upon  the  assumption  that  only  the 

relative  translational  energy  along  the  line  of  centers  at  the  time  of  impact 

is  available  for  dissociation.  Limits  can  be  placed  upon  the  value  of  Sj.  It 

can  range  in  value  from  0  when  none  of  the  internal  degrees  contribute  to 

2 

the  dissociation  energy  up  to  an  upper  limit  of  3.  Byron  found  in  oxygen 

that  Sj  was  2.  0  by  a  best  fit  of  the  data.  He  had  sufficient  sensitivity  to  show 

that  a  value  of  2.  0  for  Sj  gave  a  better  fit  of  the  experimental  data  than  a 

value  of  2.  5.  Since  the  atoms  cannot  contribute  to  the  dissociation  energy 

from  rotational  and  vibrational  degrees,  the  values  of  and  should  range 

between  the  limits  of  0  to  1.  5.  The  upper  limit  of  1.  5  arises  from  the  fact 

that  at  least  one  of  the  rotational  degrees  of  freedom  of  the  molecule  is  not 

free  to  contribute  to  the  dissociation  energy  in  order  to  conserve  angular 

momentum.  However,  an  excited  atom  can  have  a  larger  value  for  S. 

The  limitations  of  such  a  simple  theory  of  reaction  rates  have  already 

1, 2 

been  fully  discussed  in  the  literature.  It  would  be  ideal  to  have  static 
sources  for  the  study  of  high  temperature  reactions.  But  at  present  there  is 
no  such  source  available  for  temperatures  in  the  vicinity  of  5000°K.  Even  if 
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one  had  such  a  source,  so  that  very  careful  reaction  rate  measurements  could 
be  made,  this  knowledge  alone  would  not  be  sufficient  for  a  full  understanding 
of  a  shock  transition.  The  principal  motivation  of  this  research  is  toward 
a  better  understanding  of  shock  structure  in  a  diatomic  gas.  Accordingly 
then,  the  philosophy  has  been  one  of  compromise  wherein  shock  reaction  zone 
data  has  been  employed  to  infer  values  for  the  specific  dissociation  rates  in 
nitrogen.  It  is  hoped  that  the  insertion  of  these  rates  into  a  shock  structure 
model  which  includes  the  differential  equations  of  fluid  motion,  relaxation 
equations,  and  diffusion  effects  will  lead  to  more  insight  into  shock  front 
structure. 
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EXPERIMENTAL  PROCEDURE 


The  shock  tube  and  associated  equipment  used  for  this  study  are  shown 
in  Figure  1.  The  stainless  steel  tube  has  a  cross  section  of  2”  x  2".  The 
driven  part  of  the  tube  is  about  25  feet  long  while  the  driver  (not  shown)  has 
a  length  of  5  feet.  The  driver  section  is  equipped  with  5  spark  plugs  distributed 
alternately  along  the  tube  walls.  In  addition.  5  fill  ports  are  uniformly  dis¬ 
tributed  along  the  driver  for  the  purpose  of  loading  the  combustible  mixture. 

The  mixture  used  was  10%  O^.  20%  H^.  and  70%  He  at  initial  load  pressures 
of  about  300  psi,  600  psi.  and  800  psi  absolute.  The  combustible  gases  were 
allowed  to  mix  for  8  minutes  prior  to  ignition.  Water  pumped,  commercial 
grade  gases  were  used  throughout  for  the  driver  loads. 

Station  1  in  Figure  1  shows  a  light  screen  which  was  used  to  trigger  a 
tungsten  exploded  wire  light  source  (4).  If  the  light  screen  photomultiplier 
fed  directly  into  a  resistive  load,  appreciable  pre-triggering  occurred  owing 
to  the  luminosity  of  the  strong  shocks  generated.  In  order  to  reduce  this 
effect  an  inductive  load  was  inserted,  so  that  the  light  screen  was  more  sensitive 
to  high  frequencies.  Triggering  by  a  shock  roughly  corresponds  to  a  high 
frequency  signal  input  whereas  triggering  by  the  luminosity  of  the  approaching 
shock  wave  might  be  considered  to  correspond  to  a  lower  frequency  signal 
input.  In  practice,  the  inductive  load  kept  the  pre-triggering  of  the  light 
screen  to  a  maximum  of  20-30  microseconds  for  all  shock  strengths  generated 

in  this  study. 
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As  Figure  1  shows,  both  beams  of  the  Mach-Zehnder  interferometer  pass 
through  the  test  section  similar  to  Byron's^  arrangement.  The  distance  down- 
stream  from  the  light  screen  (1)  to  the  upstream  beam  of  the  interferometer 

(2)  is  2  feet  and  the  distance  downstream  to  the  other  interferometer  beam 

(3)  from  (2)  is  1  foot.  At  (2)  there  is  a  window  port  in  the  top  center  of  the 
test  section  in  which  a  Kistler  SLM  pressure  gauge  could  be  mounted. 

A  tungsten  wire  (4)  exploded  by  the  discharge  of  a  delay  line  which  had 

16  kilojoules  of  stored  energy  provided  the  light  source  for  the  recording  of 

time-resolved  inte rferograms.  This  source  has  a  rise  time  of  50  microseconds 

and  a  useful  duration  of  about  800  microseconds.  By  actual  spectrographic 

observation,  narrow  wavelength  bands  could  be  isolated  with  the  aid  of  Bausch 

o 

and  Lomb  interference  filters.  These  bands  were  never  wider  then  100  A  . 
Consequently,  good  monochromatic  fringes  could  be  obtained  in  the  blue,  green 
and  red  regions.  The  time-resolution  was  accomplished  by  use  of  a  60  micron 
slit  imaged  at  the  windows  at  locations  (2)  and  (3)  in  conjunction  with  an  air 
turbine  driven  camera.  A  beam  splitting  arrangement  inside  the  rotating  drum 
camera  permitted  the  simultaneous  recording  of  white  and  monochromatic 
fringes.  The  60  micron  slit  width  gave  good  space  resolution  and  the  ultimate 
time  resolution  was  .  15  microseconds.  All  the  interferograms  were  read  on 
an  optical  comparator  which  has  a  linear  accuracy  along  both  axes  of  1  part 
in  14.000.  Distances  could  be  measured  accurately  down  to  10  microns. 
Generally,  the  drum  was  rotated  so  that  the  film  had  a  time  scale  of  50 
microseconds  per  cm  (the  limit  is  30  microseconds  per  cm).  Ansco  70  mm 
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super,  hy-pan  film  was  employed.  Since  this  film  has  100  grains  per  millimeter, 
it  is  concluded  that  an  actual  time  resolution  of  about  .  25  microseconds  was 
attained  in  laboratory  coordinates. 

The  shock  velocity  was  measured  over  the  1  foot  distance  from  the  up¬ 
stream  beam  (2)  to  the  one  downstream  (3)  from  the  interferogram.  Simul¬ 
taneously.  the  shock  velocity  was  measured  over  the  2  feet  from  the  light 
screen  (1)  to  the  Kistler  pressure  gauge  at  (2).  The  optical  method  was 
accurate  to  .  2%  while  the  pressure  gauge  method  has  an  accuracy  ~ 

All  runs  indicated  that  the  maximum  local  attenuation  was  no  larger  than 
1%  of  a  Mach  number  per  foot  which  was  surprisingly  good.  However  this 
attenuation  agrees  well  with  that  observed  by  Wittliff  and  Waldron  for 
constant  volume  combustion  about  20  feet  downstream  of  the  diaphragm. 

Stainless  steel  and  copper  diaphragms  were  used  with  about  30%-40%  scribe 
depth.  The  x  scribes  were  held  to  a  tolerance  of  +  .001"  and  -.  000",  The 
initial  combustible  load  was  carefully  matched  to  yield  a  peak  pressure  equal 
to  the  burst  pressure  of  the  diaphragm.  A  pressure  gauge  was  mounted  in 
the  driver  for  the  purpose  of  attaining  a  good  match  between  load  and  diaphragm. 
The  material  for  the  diaphragm  was  carefully  controlled  and  the  scribe  marks 
were  cut  45°  to  the  grain.  Very  reproducible  breaks  were  obtained.  One 
series  of  3  runs  repeated  under  the  same  conditions  indicated  that  the  shock 
velocity  at  the  test  section  could  be  reproduced  to  about  2%. 

After  acquiring  good  control  of  the  constant  volume  technique,  it  was  still 

a  moot  question  whether  the  flow  region  behind  the  shock  front  would  be  steady 

2 

enough  in  time  to  permit  a  quantitative  study  of  the  reaction  zone.  Byron 
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and  Rose  have  observed  a  rapid  decay  in  gas  density  behind  the  shock  wave 

as  it  passed  a  given  point  even  when  it  was  certain  that  a  shock  and  not  a 

1 8 

detonation  wave  had  been  generated.  Byron  attributed  this  decay  to  the 

arrival  of  many  expansion  waves  from  the  hot  driver  gases  which  expand 

17 

down  the  tube.  Alternatively,  Rose  emphasized  the  importance  of  the  cooling 

of  the  shock  heated  sample  gas  by  rapid  heat  transfer  to  the  walls  through 

the  turbulent  boundary  layer.  Another  possible  source  of  these  expansion 

waves  is  the  continued  burning  of  the  driver  gases  as  they  expand  along  the 

15 

driver  tube  approaching  constant  pressure  burning  .  Some  evidence  was 
found  to  support  this  view.  On  two  runs  the  diaphragms  were  defective  and 
opened  long  before  the  burning  was  complete.  The  interferograms  seemed 
to  show  a  larger  decay.  Consequently,  it  was  believed  that  by  exercising  great 
care  in  matching  loads  and  diaphragms,  and  also  in  the  fabrication  of  the 
diaphragms,  this  decay  might  be  reduced.  The  first  interferograms  taken 
with  good  combustion  control  showed  a  promisingly  small  decay.  But  it  was 
decided  that  a  further  check  was  necessary  before  the  flow  could  be  assumed 
steady  enough  for  quantitative  studies.  A  necessary  condition  for  steady  flow 
is  that  the  reaction  zone  density  profile  be  unchanging  in  time  as  it  proceeds 
past  the  slit  at  the  upstream  interferometer  window  (2).  In  order  to  demonstrate 
that  the  flow  fulfilled  this  requirement  two  slits  placed  1  cm  apart  were  em¬ 
ployed.  The  anticipated  reaction  zone  length  prompted  this  choice  of  slit 
spacing.  In  order  to  avoid  double  photographic  exposures  the  upper  half  of  the 
upstream  slit  was  masked  off  simultaneously  with  the  lower  half  of  the  down¬ 
stream  slit.  Figure  2  shows  a  monochromatic  interferogram  taken  with  the 
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double  slit  showing  a  M  11.  0  shock  in  nitrogen  with  a  time  resolution  of  55 
microseconds  per  cm.  The  reaction  zone  can  be  clearly  seen  at  both  slits. 
Comparative  measurements  indicated  that  the  profile  shape  is  unchanged  as 
the  shock  traverses  the  1  cm  slit  separation;  hence,  it  appears  that  these 
shock  induced  flows  are  sufficiently  steady  in  time  for  measurements  to  be 
taken.  Returning  to  Figure  2,  note  that  a  vibrational  relaxation  zone  can  be 
distinguished  from  the  dissociation  zone  owing  to  the  relatively  low  Mach 
number.  A  compression  wave  can  be  seen  following  about  20  cm  behind  the 
shock.  This  wave  was  primarily  observed  only  for  combustion  runs  in 
nitrogen.  Combustion  runs  into  the  mixtures  and  cold  hydrogen  pressure 
breaks  into  nitrogen  did  not  generate  any  sizeable  secondary  waves. 

The  presence  of  such  waves  pursuing  the  shock  wave  along  the  tube 
attests  to  the  desirability  of  examining  the  flow  field  behind  combustion 
driven  shock  waves  prior  to  making  other  types  of  measurements,  for  example, 
spectrographic. 

The  fringe  shifts  could  be  read  accurately  to  at  least  .  1  of  a  fringe  width. 
However,  in  adjacent  regions  along  the  profile  the  fringes  were  sensitive  to 

displacements  of  only  .03  fringe  which  agrees  with  the  conclusion  reached  by 

2 

Mathews  . 

Figure  3  shows  a  M  10.3  wave  in  a  13%  nitrogen-87%  neon  mixture.  For 
simplicity  only  the  white  interferogram  is  shown.  When  the  shock  arrives  at 
the  upstream  window  the  fringes  jump  up  corresponding  to  the  rise  in  density. 
Then  there  is  a  small  decay  which  is  quite  linear  until  the  shock  arrives  at  the 
downstream  window.  Here  the  fringes  jump  back  below  the  original  position, 
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because  the  density  jump  turns  out  to  be  exactly  the  same  as  it  was  when  the 
shock  was  at  the  upstream  window.  This  is  consistent  with  the  shock  velocity 
measurements  which  indicated  that  there  was  no  more  than  a  1%  decrease  in 
Mach  number  of  the  shock  as  it  traversed  the  one  foot  distance  between  the 
two  windows.  The  vibrational  relaxation  zone  can  also  be  seen  in  Figure  3. 

Figure  4  shows  a  monochromatic  interferogram  of  a  M  11.5  shock  in  a 
15%  nitrogen-85%  argon  mixture.  The  decay  is  again  relatively  small.  Through¬ 
out  this  study  the  decay  was  generally  small  enough,  so  that  its  effect  was 
negligible  on  the  reaction  zone.  However,  in  the  neon  mixtures  the  decay 
began  to  have  some  influence  as  equilibrium  was  approached.  Therefore  in 
these  runs  the  reaction  profile  could  only  be  followed  up  to  about  70%  of 
equilibrium.  The  interferometric  sensitivity  was  smallest  in  the  neon  mix¬ 
ture.  To  partially  offset  this  fact,  all  the  runs  with  the  neon  mixture  were 
made  at  initial  pressures  of  about  2  cm. 

Fundamentally,  the  interferometer  measures  changes  in  the  index  of 
refraction  of  the  gas.  These  can  be  linearly  related  to  density  changes  through 
the  Gladstone- Dale  relation  provided  the  mean  specific  refractivity  is  known. 
White has  measured  the  specific  refractivity  of  the  nitrogen  atom  and  found 

19  . 

good  agreement  with  a  calculation  made  by  Dalgarno  .  Since  he  used  an  in¬ 
terferometer  it  may  be  argued  that  perhaps  dissociative  equilibrium  is  not 
obtained  behind  the  shock.  If  Dalgarno' s  calculation  is  good  to  an  order  of 
magnitude,  and  it  appears  to  be  so,  then  one  can  select  degrees  of  equilibrium 
dissociation  for  which  the  specific  refractivity  of  the  nitrogen  atom  will  play 
a  negligible  role  as  long  as  its  order  of  magnitude  is  known.  For  example, 
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if  the  shock  leads  to  only  6%  equilibrium  dissociation,  then  the  index  of 

refraction  measurements  can  be  directly  converted  into  densities  safely 

neglecting  the  specific  refractivity  of  the  nitrogen  atom.  In  this  way,  ar 

unambiguous  demonstration  can  be  made  showing  that  equilibrium  is  attained 

at  lower  Mach  numbers  (and  therefore  temperatures);  it  is  even  more  certainly 

attained  at  higher  Mach  niunbers  where  for  higher  degrees  of  dissociation 

(10%  and  upwards)  one  must  take  into  account  the  nitrogen  atom  refractivity. 

3 

White's  measured  value  of  .  31  —  for  the  specific  refractivity  of  the 

gm 

nitrogen  atom  was  incorporated  into  the  theoretical  profiles  of  the  index  of 
refraction  as  a  function  of  distance  behind  the  shock  front  in  shock  fixed 
coordinates.  The  distances  measured  on  the  film  were  converted  to  distances 
in  shock  fixed  coordinates  by  the  following  relation 

(shock  fixed)  = 

where  U  is  the  shock  velocity  and  x  is  the  distance  measured  along  the  time 
s 

axis  of  the  film  from  the  shock  front.  The  film  speed  which  was  measured  by 
a  Berkeley  counter  and  the  shock  velocity  were  known  to  an  accuracy  of  -j—  %. 
This  conversion  introduced  no  error  since  owing  to  the  limitations  of  film 
grain,  measurements  taken  in  shock  fixed  coordinates  would  not  be  better 
than  %.  Film  shrinkage  after  development  was  monitored  and  found  to 
introduce  negligible  error  in  the  measurements. 

It  will  be  instructive  to  consider  the  actual  magnitude  of  the  numbers 
involved  in  two  cases  in  order  to  gain  insight  into  the  influence  of  changing 
refractivity  versus  changing  density  upon  the  index  of  refraction. 
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Case  I:  100%  nitrogen 

Mach  No.  =  13.  5;  9.95%  of  the  molecules  dissociated  at  equilibrium 
The  Gladstone-Dale  relation  is 

(n  -  1)  =  bp  (5) 

where  b  is  the  specific  refractivity,  p  the  density,  and  n  the  total  index  for  any 

point  on  the  reaction  profile.  Since  the  interferometer  measured  changes  in 

total  index  let  b  ,  p  ,  and  n  represent  the  same  quantities  at  the  point 
o  o  o 

where  vibration  is  equilibrated  and  no  dissociation  has  begun.  The  quantities 

in  Eq.  (5)  may  be  defined  as 

n  =  n  +  An 
o 

p=p+Ap 

o 

b  =  b  +  Ab 

o 

Substituting  into  (5) 


(n  -  1)  +  An  =  b  p  + 
0  00 

p  Ab  +  b  A  A  +  Ao  Ab 

0  0  ^ 

.  *  .  An  =  b  Ap  +  Ab 

0 

[p„  ♦  Ap] 

(7) 

Since  n  -l  =  b  p  by  (5) 
o  o  o 

The  change  An  from  the  point  corresponding  to  vibrational  equilibrium  and  no 

dissociation  up  to  full  equilibrium  is  8.  14  x  10  .  The  corresponding  values 

of  b  ,  Ab,  An,  and  n  are, 
o  ^  *^0 

3  3 

(  X  =  4549  A)  :  b  =  .  241  ;  Ab  =  6.50xl0"  ;  ;  =  30.28x10  gm/cm 

'  o  gm  gm 

p  =  98.  51X10  ^  gm/cm^ 
o 

Substituting  into  Eq.  (7) 
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8.  14  X  10"^  =  (.  241)  (30.  28)  lO"^  +  (6.  50)  lO"^  [98.  51  x  10 


+  30.  28  X  lO"^]  =  7.  30  X  lO"^  +  .  84  x  lO"^ 

Note  that  the  b  Ap  term  is  about  9  times  larger  than  the  term  involving  Ab. 
o 

Case  ir.  15%  nitrogen-85%  argon 

Mach  No.  =  11.  7;  80%  of  the  molecules  dissociated  at  equilibrium 

The  differences  between  the  values  at  full  equilibrium  and  vibrational 
equilibrium  are: 

3  cm^ 

An  =  15.  10  X  10"^;  Ap  =  77.  05  x  10‘  gm/cm  ;  Ab  =  6.  26  x  10 
Substituting  into  Eq.  (7): 

15.  10  X  lO”^  =  (.  169)  (77.  05)  10‘^  +  (6.  26)  10‘^  [252.  77xl0"^ 

+  77.05  X  10"^]  =  13.03  x  lO"^  +  2.  07  x  lO'^ 

Consequently,  while  the  first  term  on  the  right  still  dominates  it  is  now  only 
about  6  times  larger  than  the  second  term  on  the  right.  If  the  error  in  White's 
value  were  as  large  as  5%  it  would  still  only  produce  a  1%  error  in  the  change 
in  total  index  in  Case  II  above. 

The  initial  pressure  of  the  sample  gas  prior  to  shock  heating  was 
measured  by  a  Wallace  and  Tiernan  absolute  manometer,  and  pressures  above 
2  cm  of  mercury  were  read  on  an  absolute  mercury  manometer  with  a 
cathetometer.  The  manometer  was  isolated  from  the  system  by  a  liquid 
nitrogen  trap.  The  initial  temperature  of  the  gas  was  measured  with  a 
calibrated  thermometer  which  could  be  read  directly  to  tenths  of  a  degree. 
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In  regard  to  the  pre-mixing  of  driver  and  test  gas  by  diffusion  across  the 
contact  zone,  Hooker  has  shown  by  spectrographic  absorption  measure¬ 
ments  for  initial  gas  pressures  of  .  5  cm  and  upwards  that  40%  of  the  ideal 
test  time  is  free  from  driver  gas  contamination.  The  initial  gas  pressure 
ranged  from  about  9  mm  -  27  mm  in  the  present  study.  In  all  of  the  interfero- 
grams  a  sharp  rise  in  index  of  refraction  could  be  seen  entering  the  upstream 
beam  of  the  interferometer  after  the  shock  had  already  arrived  at  the  down¬ 
stream  beam.  In  Figures  3  and  4  this  is  easily  seen  and  probably  represents 
the  arrival  of  the  driver  gases.  In  all  cases  the  actual  test  time  appeared  to 
be  from  33%-50%  of  the  ideal.  The  Kistler  pressure  gauge  records  indicated 
that  the  average  value  of  the  pressure  was  constant  over  this  whole  zone  of 
flow  up  to  50%  of  the  ideal  flow  duration.  As  far  as  boundary  layer  effects 
are  concerned  they  are  negligible  in  all  cases  over  the  length  of  reaction  zones 
observed  here  (1  cm-3  cm).  Interferograms  of  shocks  in  pure  argon  at  Mach 
numbers  of  about  6  conveniently  reveal  information  on  rate  of  growth  of  the 
boundary  layer. 

Throughout  this  study  Matheson  pre-purified  water  pumped  nitrogen  and 
argon  were  used.  These  were  guaranteed  to  have  impurity  levels  less  than 
20  ppm.  The  neon  was  of  research  grade.  In  the  early  pure  nitrogen  runs 
the  system  was  evacuated  to  about  1  micron  and  in  the  later  runs  to  .  6  micron. 
This  was  followed  by  a  continuous  flush  at  the  initial  pressure  chosen  for  each 
run  up  to  within  20  seconds  of  ignition.  The  total  leak  rate  averaged  between 
1-2  microns  per  minute.  When  the  mixtures  were  run  a  larger  diffusion  pump 
was  put  into  the  system  and  the  system  was  pumped  down  to  about  .  2  micron 
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read  on  a  Phillips  gauge.  All  diffusion  pumps  were  isolated  by  liquid  nitrogen 
traps.  To  further  check  the  cleanliness  of  the  system,  a  glass  tonometer 
tube  was  fuzed  into  the  tube.  A  typical  pre-run  cycle  was  carried  out  and  then 
the  tonometer  vessel  was  brought  up  to  atmospheric  pressure  with  the  sample 
gas  and  sealed  off.  Crippen  laboratories  made  a  mass  spectrometric  analysis. 
They  found  that  oxygen,  hydrocarbons,  and  water  vapor  were  present  in 
quantities  less  than  25  ppm  which  is  in  good  agreement  with  levels  expected 
from  a  knowledge  of  the  vacuum  system  and  gas  purity. 

A  spectr ographic  window  was  inserted  in  the  top  of  the  tube  above  the 
upstream  interferometer  beam  in  place  of  the  pressure  gauge.  An  f-3  time- 
resolved  spectrograph  (See  5  in  Figure  1)  was  utilized  to  examine  the  luminosity 
emitted  behind  shocks  in  nitrogen  and  the  nitrogen-argon  mixture.  The  violet 
band  system  of  CN  was  not  observed  since  it  would  have  to  be  present  at 
levels  of  about  150  ppm  in  order  to  be  observed  with  the  present  spectrograph. 
The  first  negative  band  system  of  was  not  observed  in  pure  nitrogen 
shocks  up  to  Mach  numbers  of  14. 

Since  the  combustion  technique  produces  large  amounts  of  water  in  the  tube 
it  was  felt  advisable  to  manually  clean  the  tube  prior  to  repumping.  Wooden 
blocks  covered  with  a  relatively  lint-free  paper  were  pushed  through  the  tube. 
After  this  was  done  a  careful  inspection  was  made  by  shining  a  light  down  the 
driven  section  of  the  tube  to  be  certain  that  no  dust  or  lint  remained. 
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RESULTS 


Tables  I  and  II  give  the  pertinent  data  for  fourteen  final  data  runs  yielding 

information  on  vibrational  relaxation  times  and  dissociation  rates.  The 

vibrational  relaxation  times  obtainea  in  pure  nitrogen  were  about  a  factor 

20 

of  two  longer  than  Blackman's  .  Previous  experimentation  has  shown  that 

the  vibrational  relaxation  time  can  be  shortened  considerably  if  there  are 

small  amounts  of  water  vapor  and  hydrocarbons  present.  The  longer  times 

observed  here  may  be  a  direct  result  of  better  impurity  control.  The  pure 

nitrogen  data  of  Table  I  do  not  cover  a  sufficient  temperature  range  to 

21  20 

distinguish  between  Widom's  and  Landau's  theories. 

Turning  from  the  results  for  pure  nitrogen  to  the  mixtures,  the  neon 
mixture  dissociation  runs  had  a  surprisingly  long  vibrational  relaxation  zone. 
Within  the  precision  of  the  measurements  neon  seemed  to  be  less  than  10%  as 
effective  as  nitrogen  molecules  in  exciting  nitrogen  vibrationally.  The  argon 
mixture  data  given  in  Table  I  implied  that  argon  was  about  as  efficient  as 
the  nitrogen  molecule  in  exciting  nitrogen  vibrationally.  The  argon  efficiency 
was  inferred  by  comparing  the  actual  relaxation  time  in  the  mixture  to  an  extra¬ 
polated  time  based  on  the  lower  temperature  nitrogen  data, assuming  that  the 
argon  played  no  role.  Extrapolation  by  either  Widom's  or  Landau's  theory 
inferred  that  the  argon  atom  is  as  efficient  in  vibrationally  exciting  nitrogen 
as  the  nitrogen  molecule  itself.  This  result  checks  well  with  the  observation 

by  Blackman  that  the  nitrogen  molecule  was  only  40%  as  effective  as  the  oxygen 

18 

molecule  in  exciting  oxygen  vibrationally.  Byron  found  that  the  argon 
atom  was  about  30%  as  efficient  as  the  oxygen  molecule  in  exciting  oxygen 
vibrationally. 
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Concerning  the  question  of  vibrational  equilibration  in  the  dissociation 
runs  the  following  procedure  was  adopted.  First,  the  time  was  measured 
behind  the  shock  front  to  the  point  which  had  the  index  of  refraction  correspond¬ 
ing  to  vibrational  equilibrium.  This  time  was  then  compared  with  the  relaxation 
time  extrapolated  by  the  Landau  theory  (predicting  slightly  longer  times  that 
Widom's),  and  in  every  case  the  observed  time  was  no  longer  than  Landau's. 

For  the  higher  temperature  runs  in  the  argon -nitrogen  mixture  the  vibrational 
relaxation  zone  was  so  short  that  it  coalesced  with  the  shock  front.  In  closing 
this  discussion  of  the  vibrational  relaxation,  it  may  safely  be  concluded  that 
vibration  was  equilibrated  before  any  measurable  dissociation  occurred. 
Therefore,  there  was  no  detectable  coupling  between  the  vibrational  and 
dissociative  processes  over  the  temperature  and  pressure  ranges  of  these 
experiments. 

Referring  to  Table  II  it  is  seen  that  the  highest  temperature  obtained  for 
vibrational  equilibrium  and  no  dissociation  was  10,  963°K  (no  reaction).  The 
experimentally  inferred  rates  imply  that  even  in  this  case  only  about  1  in 
10,000  collisions  produces  a  dissociation.  Therefore,  the  assumption  that 
the  vibration  and  rotational  degrees  maintain  equilibrium  with  the  local 
kinetic  temperature  throughout  the  dissociation  process  is  valid,  since  it 
takes  only  about  1,000  collisions  to  equilibrate  vibration,  and  less  than  100 
collisions  to  equilibrate  rotation  at  these  temperatures.  Consequently,  use  of 
the  collision  rate  expressions  given  in  Eq.  (3)  which  ignore  coupling  between 
the  dissociation  and  vibrational  processes  is  justified  in  this  study. 
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The  initial  values  that  were  chosen  for  S^,  S^.  and  in  Eq.  (2)  were  2.0, 
1.0,  and  1.0  respectively.  These  were  the  values  reported  by  Byron^  from 
data  taken  in  nitrogen-argon  mixtures.  The  IBM  1620  could  run  twelve 
different  profiles  simultaneously  for  various  choices  of  ♦j,  ^  ,  and  A 

good  over-all  fit  of  the  nine  runs  of  Table  I  was  not  obtained  for  these  values 
of  S.  Figures  5,  6,  7,  and  8  show  the  experimental  profiles  for  runs  3,  5,  6, 
and  7.  A  theoretical  profile  is  shown  in  each  of  the  figures  which  has  the  S 
and  ♦  values  obtained  by  Byron.  A  glance  at  Table  II  shows  that  runs  3,  5,  6, 
and  7  cover  almost  a  3,  000°K  range  in  the  no  reaction  temperature.  Further, 
owing  to  the  higher  dissociation  energy  of  nitrogen,  the  temperature  drop 
between  the  no  reaction  point  and  dissociative  equilibrium  is  greater  than  that 
occurring  in  oxygen  for  the  same  degree  of  equilibrium  dissociation.  In  view 
of  the  fact  that  a  good  fit  could  not  be  obtained  for  Byron  s  values  of  Sj,  S^. 
and  Sj,  it  was  decided  to  vary  these  parameters.  The  value  of  Sj  was  varied 
between  1.7  and  2.5  in  incremental  steps.  The  corresponding  values  of 
and  were  set  equal  to(Sj  —  1.0)  in  each  case. 

Then  twelve  profiles  were  generated  for  runs  3,  6,  and  7  for  each  set  of 
S  values  by  varying  the  f  values.  An  optimum  set  of  f  values  could  be  picked 
for  every  set  of  S  values.  After  plotting  the  profiles,  the  following  set  appeared 
to  hold  promise  for  obtaining  a  good  compromise  fit  for  the  profiles  of 
Figures  5,  6,  7,  and  8; 

I 

S,  =  2.2;  S_  =  S,  =  1.2;  #  ,  =  .08  ;  =  .3;  =  .08 

12  3  c  j 

At  this  point  it  was  discovered  that  there  was  an  erroneous  factor  of  V3~ 
multiplying  the  expression  for  k^^.  Therefore  the  correct  value  of 
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that  cited  above  multiplied  by  1.73;  hence,  in  these  three  cases  the  true  values 
of  0^  larger  than  the  values.  This  was  a  disturbing  result,  since  if 

the  classical  collision  theory  has  any  validity  at  all,  it  is  expected  that  the 
efficiency  parameter  for  collisions  with  rare  gas  atoms  should  be  no  larger 
than  which  pertains  to  collisions  with  nitrogen  molecules.  Therefore  it 
was  decided  to  vary  and  (letting  =  S^)  while  was  held  fixed  at  values 
of  1.7,  2.0,  2.2,  2.3,  2.4,  and  2.  5.  The  argon-nitrogen  mixtures  offered  the 
best  chance  for  detecting  any  sensitivity  to  variation  of  and  since  they 
had  the  largest  degrees  of  dissociation  and  covered  a  temperature  range  of 


3.  000°K.  With  the  values  of  0.'  =  .08  and  S  =  1.3  the  magnitude  of  the 


product  4 


f-Y 

i  \RT/  * 


^V3~was  determined  for  the  hottest  run  (6).  Recalling  that 


=  ♦3'  VT.  as  eries  of  ^  and  S  values  were  calculated  that  kept  the  value 


of 


i(^)  ’ 


so  that  a  good  fit  of  the  upper  part  of  the  profile  in  Figure 


7 (run  6) was  maintained.  These  choices  were  then  fitted  to  the  coolest  nitrogen- 


argon  run  in  Figure  6;  The  profile  seemed  to  favor  the  following. 


Sj  =  2.2;  =  1.5;  =  .08  ;  ^3  =  .08 

This  fits  slightly  better  at  the  top  of  the  profile  than  does  Sj  =  2.2;  =  =  1.2  ; 

^  =  .08  ;  ♦  =  .3  :  ^  =  .  18.  Admittedly,  there  is  no  great  precision  in 
1  2  3 

getting  the  best  value  for  and  S3;  However,  the  data  suggested  that 
the  values  of  a^nd  S3  were  nearer  to  1.  5  than  1.  0.  The  value  adopted 
here  for  S^  and  S3  of  1.  5  is  still  consistent  with  the  predictions  of  classical 
collision  theory. 
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In  regard  to  the  determination  of  S^,  Figures  5,  7,  and  8  show  that  2.3 
was  slightly  better  than  2.  0  or  2.  5.  The  final  set  of  values  adopted  here 
was: 


Sj  =  2.2;  =  =  1.5;  =  .08  =‘^5;  =  .08 

Hence  the  final  values  for  the  rate  constants  inferred  here  are: 


(N,  -N,):k^  =  (6.8x10^^  X  x  exp'(^^V^°  )±  25% 


(N,  -  N  )  :  =  (5.5x  10^*^)  X  t"^‘°  x  exp  ( 

2 


T 

-,113,  260 


mole-sec. 


)±  40% 


cc 


mole-sec. 


.19.  ..,-1.0 


(N.,  -  Ne)  ;  k^  =(1.5x10  )xT 
2  d3 


X  exp  (• 


113,260 


)  ±  40% 


cc 


mole-sec. 


(N  -A)  )  :  k  =  (2. 1  X  10^^)  X  T  '''  .  x  exp  (  ±  40% 

4 


mole-sec. 


The  ratio  of  the  two  values  for  neon  and  argon  as  third  bodies  was 
determined  from  the  conservation  equations  of  energy  and  momentum.  The 
fraction  of  energy  transferred  to  a  nitrogen  molecule  on  impact  was  calculated 
for  both  rare  atoms.  Equating  the  initial  kinetic  energies  of  the  two  rare 
atoms  before  they  undergo  impact  with  the  nitrogen  molecule  allows  one  to 
estimate  a  relative  efficiency  for  each  in  transferring  energy  to  the  nitrogen 
molecule.  Argon  and  neon  turned  out  to  be  equally  efficient  by  this  calcula¬ 
tion  and  so  was  taken  to  have  the  same  value  for  both  argon  and  neon. 

This  estimate  of  efficiencies  is  purely  classical  and  rests  upon  the  assumption 
that  the  rare  atoms  can  be  treated  as  rigid  spheres.  This  is  a  fair  assumption 
for  neon  owing  to  its  high  excitation  potentials,  but  was  somewhat  weaker  for 
argon.  An  examination  of  the  luminosity  behind  the  shocks  in  the  nitrogen- 
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argon  mixtures  by  means  of  the  f-3  spectrograph  indicated  no  evidence  of 
argon  excitation. 

Table  III  gives  the  values  of  the  inferred  recombination  rates  at  6,  400  K 

for  the  four  third  bodies  (N^.  N.  A.  Ne).  Byron's^  values  with  N^.  N.  and 

4 

A  as  third  bodies  are  shown  also,  as  well  as  Eckerman's  value  with  Xe  as 

the  third  body,  and  Camm's^^  value  for  and  N  as  third  bodies.  Recently 

Allen^  has  collected  more  data  using  the  same  technique  as  Camm.  While 

the  value  for  the  recombination  rate  constant  k  is  still  about  the  same  with 

N  as  the  third  body  the  value  of  k  is  at  least  13  times  smaller  than  k^  . 

1  2 
This  result  does  not  seem  reasonable.  The  success  of  their  method  rests 

upon  the  validity  of  the  kinetic  model  for  the  origin  of  the  radiative  decay 

that  is  observed.  Hammerling's  proposed  kinetic  scheme  has  been 

criticized^^.  This  is  an  ingenious  experiment,  but  it  appears  that  there  is 

need  for  further  experimental  studies  of  the  origin  of  the  luminosity  behind 

strong  shocks  in  nitrogen. 

In  addition  to  presenting  the  inferred  rate  constants  at  6,400  K  in 
Table  HI,  the  extrapolated  values  for  k^  (N^  third  body)  and  k^  (A  third 
body)  adopted  in  this  paper  and  those  of  Byron  are  listed  for  comparison  with 
the  measurements  of  Harteck  ,  Heron  ,  and  Wentick  made  at  300  K.  The 
rate  constant  k,  obtained  in  this  paper  is  about  13-20  times  larger  at  300  K 
than  the  other  values.  This  failure  of  rate  constants  inferred  from  shock  tube 
measurements  to  extrapolate  down  to  the  values  obtained  at  300  K  has  caused 
criticism  of  the  shock  tube  measurements.  Therefore,  some  discussion  of 
this  point  is  warranted.  If  we  wished  to  choose  the  appropriate  value  of 
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in  the  expression  for  k,  at  300°K,  we  must  consider  how  many  internal  degrees 
of  freedom  can  contribute  to  the  dissociation  of  one  nitrogen  molecule  upon 
with  another*  Since  the  vibrational  characteristic  temperature  of 
nitrogen  is  high  (3,337°K),  only  about  1  in  66,000  nitrogen  molecules  will  possess 
the  full  classical  vibration. 

Classically,  probably  only  2  rotational  degrees  and  the  vibrational  degree 
of  the  dissociating  molecule  can  contribute  to  the  energy  of  dissociation;  hence, 

Sj  takes  on  a  value  of  2.  0.  Then  by  making  a  reduction  in  0^,  from  .  08  to  .  03 
and  setting  equal  to  2.  0,  the  calculated  rate  k^^  will  agree  fairly  well  with  the 
measurements  of  Herron,  Wentinck,  and  Harteck.  Since  the  collision  theory 
is  only  classical,  it  is  unreasonable  to  demand  that  and  S^  remain  constant 
over  a  temperature  range  of  6,  000°K.  The  changes  in  S^  and  suggested 
here  (.  2  and  .05)  be^veen  6,  400°K  and  300°K  are  not  totally  inconceivable.  At 
the  high  temperatures,  an  S^  of  2.2  seems  to  give  a  better  fit  than  an  S^  of 
2.  0.  If  the  vibrationil  energy  of  both  colliding  molecules  were  available  for 
dissociation,  S^  would  be  2.5.  The  fact  that  a  value  of  2.2  seems  to  be 
favored  at  the  high  temperatures  is  therefore  not  unreasonable,  implying 
that  some  of  the  vibrational  energy  of  the  undissociating  molecule  is  available. 
Recall  that  over  the  temperature  range  of  these  experiments  almost  all  of  the 
nitrogen  molecules  possess  their  full  classical  vibration. 

Another  criticism  of  shock  tube  dissociation  measurements  is  the 
reported  anonomously  high  efficiency  for  dissociation  of  the  parent  molecules 
by  the  daughter  atoms.  All  of  the  shock  tube  measurements  in  oxygen  and  the 

24 

recent  work  of  Gardiner  in  hydrogen  have  yielded  this  result.  The  present 
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study  differs  from  those  measurements  in  that  the  results  here  imply  that  the 
nitrogen  atom  is  about  as  efficient  as  the  nitrogen  molecule  in  promoting  the 
dissociation  of  nitrogen  molecules  at  6,400  K,  It  should  be  emphasized  that  the 
different  nitrogen  atom  rate  constant  found  here  is  not  a  direct  consequence 
of  the  fact  that  the  data  favors  values  nearer  to  1.  5  for  S_  and  S,  than  1.0, 

/  \  g  ^ 

The  rate  is  proportional  to  the  product  *  Therefore,  the  fact 

that  k  is  of  the  same  order  as  k ,  at  6,  400  K  is  independent  of  the  ability 
•^2  '*1 

to  select  S-  or  S  precisely.  Ihe  fact  that  there  seemed  to  be  a  tendency  of  the 
data  to  favor  nearer  to  1.  5  than  1.  0  can  be  attributed  to  the  use  of  the 
combustion  technique  which  allows  one  to  achieve  high  degrees  of  dissociation 
(80%)  at  high  initial  pressure  (3  cm)  for  mixtures  that  are  not  too  lean  (15% 
N2-85%  A).  This  fact  helps  the  sensitivity  of  any  optical  method  of  study. 

A  comparison  of  the  sensitivity  of  the  measuring  techniques,  of  tube  widths, 
initial  pressures,  composition,  and  percent  change  in  density  throughout  the 
reaction  zone  for  the  same  degree  of  equilibrium  dissociation  between  this 
study  and  that  of  Rink  indicates  that  there  was  3-4  times  the  over-all 
sensitivity  in  this  study.  It  is  not  to  be  implied  that  the  results  here  cannot 
be  improved  upon.  More  data  and  higher  Mach  number  shocks,  especially 
in  the  neon  mixtures  at  higher  initial  pressures  of  about  3  cm,  would  be  of 
considerable  aid  in  determining  more  precisely  values  for  S^,  S^,  S^,  and 


^1'  ^2'  ^3- 

One  final  remark  should  be  made  in  connection  with  the  curve  fitting 

procedure.  The  neon  mixture  runs  were  more  sensitive  to  variations  in 

and  ^  ,  S  than  ^  ,  S  .  Whereas  in  the  argon  mixture  the  initial  part  of  the 
3  3  2  2 
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profile  was  sensitive  to  and  but  not  too  sensitive  to  or  the  upper 
part  of  the  curve  was  very  sensitive  to  and  The  pure  nitrogen  runs 
were  very  sensitive  to  and  throughout  the  entire  profile.  They  showed 

sensitivity  to  ^2  and  only  at  the  top. 

Returning  to  Table  III.  Eckerman's^  implied  high  third  body  efficiency  for 

xenon  may  possibly  be  a  direct  consequence  of  excitation  of  the  xenon,  since  its 

10 

first  excited  state  lies  at  about  8.  3  e.  v.  The  result  of  Harteck  that  the 
argon  is  equally  efficient  is  the  nitrogen  molecule  as  a  third  body  is  in 
disagreement  with  Herron's  work  and  with  Britton's  data  for  the  relative 
efficiencies  of  atoms  and  molecules  in  promoting  the  dissociation  of  iodine. 

Returning  to  Byron's  studies  in  nitrogen-argon  mixtures,  a  few  specific 
differences  should  be  pointed  out  between  his  study  and  the  present  one. 

Byron  only  matched  one  point  of  the  profile,  neglected  to  take  into  account  the 
change  in  refractivity  with  increasing  nitrogen  atom  population,  and  took  his 
measurements  from  the  central  white  fringe  shift.  Profiles  were  generated 
for  his  highest  Mach  number  run  generated  in  his  leanest  mixture  (5%  N2-95%  A) 
The  results  showed  that  his  neglect  of  the  recombination  rate  up  to  the  point  of 
apparent  50%  dissociation  led  to  no  significant  error.  Further,  calculations 
based  upon  room  temperature  dispersion  data  for  nitrogen  indicated  that  the 
overshifting  of  the  central  white  fringe  did  not  affect  the  accuracy  of  his 
measurement.  But  the  neglect  of  the  increasing  refractivity  indicated  that 
his  apparent  50%  dissociation  time  might  be  20%  too  long,  thus  causing  an 
underestimate  of  the  rate  for  his  hottest  mixture.  This  error  would  be 
negligible  for  his  cooler  runs,  so  that  it  would  not  affect  his  over-all  choice  of 
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0values.  On  the  contrary,  if  one  neglected  the  change  in  refractivity  or  read 
the  white  fringe  shifts  rather  than  the  monochromatic  shifts  for  the  argon- 
nitrogen  mixtures  in  this  study,  sizeable  errors  would  be  introduced  even  into 
the  measurement  of  the  50%  dissociation  time,  and  the  apparent  profile  would 
be  distorted.  Therefore,  the  prime  reason  for  the  results  of  this  paper  differing 
from  those  of  Byron  is  the  greater  sensitivity  here  and  the  greater  temperature 
span  covered. 

The  very  fact  that  a  theoretical  model  which  neglects  diffusion  effects 
gives  such  a  fair  fit  of  the  observed  data  over  a  wide  temperature  range  is  in 
itself  strong  evidence  that  diffusion  effects  played  no  significant  role.  For 
example,  the  principal  transport  phenomenon  is  expected  to  be  heat  conduction 
rather  than  mass  diffusion.  Heat  conduction  is  largest  at  the  steep  beginning 
of  the  profile.  Calculation  based  upon  the  maximum  initial  temperature 
gradient  at  the  beginning  of  the  reaction  profile  of  run  number  3  in  pure 
nitrogen  indicated  that  heat  transport  affected  the  index  of  refraction  profile 
by  less  that  1/2%  even  at  the  steepest  part.  Over  the  major  portion  of  the 
profile  the  effect  was  negligible.  In  the  neon  mixtures  the  gradients  were 
much  smaller,  so  that  heat  transport  played  no  measurable  role. 

Since  the  theory  assumes  that  vibration  is  equilibrated  prior  to  dissociation 
in  those  runs  where  a  vibration  zone  was  detected,  the  origin  of  the  profiles 
was  shifted  along  the  X  axis  to  the  point  at  which  vibrational  equilbrium  was 
attained.  In  Figures  6  and  7  in  the  argon-nitrogen  runs,  the  vibrational 
zone  was  coalesced  with  the  shock  front,  so  that  the  origin  was  not  shifted. 

The  relative  and  absolute  errors  in  measurement  are  indicated  on  Figures  5 
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through  8.  The  error  in  fixing  the  position  of  the  shock  front  is  shown  in  Figures 
6  and  7.  Although  only  nine  runs  were  used  to  infer  the  dissociation  rate 
measurement,  many  more  runs  were  made  in  the  over-all  study.  At  least 
a  half-dozen  were  made  using  pressure  gauges  in  the  window  ports  for  the 
interferometer  beams  in  order  to  measure  attenuation  of  the  shock.  Another 
half  dozen  interferometer  runs  were  made.  These  for  the  most  part  attested 
to  the  good  control  of  the  combustion  technique.  They  were  not  used  for  data 
reduction  for  other  reasons:  small  fringe  shifts  in  the  reaction  zone,  uncertain 
initial  pressure  owing  to  failure  of  the  vacuum  valve  to  fully  close,  underexposure 
of  the  optical  image,  and  splitting  of  the  optical  image  across  the  butted  ends 
of  the  film.  Out  of  the  twenty-six  interferometric  runs  made  with  combustion 
driven  shocks,  on  one  occasion  two  shocks  were  observed  and  the  diaphragm 
showed  signs  of  a  bad  opening.  On  another  occasion  a  wave  was  seen  overtaking 
the  shock  front  in  the  reaction  zone  itself;  however,  no  visible  evidence  was 
found  of  faulty  diaphragm  opening. 
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CONCLUSIONS 


1.  It  has  been  experimentally  demonstrated  that  combustion  generated 
shocks  may  be  used  for  shock  structure  studies  providing  that  the  peak 
burning  pressure  of  the  combustible  loads  is  carefully  matched  to  the  burst 
pressure  of  the  diaphragms.  The  flow  field  behind  the  shock  should  be  moni¬ 
tored  simultaneously  with  the  recording  of  other  types  of  data  for  the  most 
reliable  results. 

2.  Owing  to  the  greater  sensitivity  attained  by  operating  with  not  too  lean 
mixtures,  higher  initial  pressures  (1  cm-  3  cm)  and  high  shock  Mach  numbers, 
a  more  reasonable  value  was  obtained  for  the  ratio  of  the  efficiency  of  the 
nitrogen  atom  to  the  nitrogen  molecule  in  promoting  the  molecular  dissociation. 
This  ratio  of  about  unity  at  6,  400°K  is  more  consistent  with  the  results  of 
careful  room  temperature  measurements. 

3.  In  view  of  the  over-all  good  fit  of  the  data,  it  is  concluded  that  the  naive 
picture  of  vibrational  relaxation  being  completed  ahead  of  the  onset  of  disso¬ 
ciation  appears  to  be  valid  for  the  temperature  and  density  ranges  covered 
in  this  study.  Consequently,  it  would  appear  that  the  inferred  rate  constants 
here  approximate  the  true  values  quite  well  over  the  temperature  range  in 
this  study.  However  it  would  be  a  mistake  to  use  these  inferred  rate  constants 
at  higher  temperatures  where  coupling  between  the  vibrational  relaxation  and 
dissociation  cannot  be  ignored.  If  the  inferred  rate  constants  are  used  at 
temperatures  intermediate  between  room  temperature  and  the  range  here 

(5,  795°K  -  10,  963°K),  they  will  be  in  error  (not  owing  to  coupling  between 
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vibration  and  dissociation,  )  because  the  correct  values  of  and  will 
probably  be  smaller. than  those  reported  here. 
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290.9  17.57  4.  14  35.0  6515  7845 


TABLE  III 


k  cm 

Inferred  Recombination  Rates  - ^ - 

mole  -  sec 


Source 

Third  Body 

300°K 

6,  400°K 

This  paper 

^2 

1.6  .  10^^ 

4.2  • 

10^^ 

N 

2.0  .  10^^ 

4.7  • 

lo'^ 

Ne 

5.4  .  10^^ 

1.  3  • 

10^^ 

A 

8.0  .  10^^ 

1.  9  • 

lo'" 

5 

.  16 

.^14 

Byron 

'^2 

8.6  •  10 

4.  4  • 

10 

N 

1.9  •  10^^ 

2.  0  • 

10^5 

A 

1.5  •  10^^ 

1.  6  • 

10^^ 

4 

Eckerman 

Xe 

5.  1  •  10^^ 

2.  7  • 

10^5 

22 

15 

Camm 

'^2 

<1.  1  • 

10 

15 

N 

8.7  • 

10 

Allen^ 

N 

6.  5  • 

10^5 

8 

Herron 

^2 

5. 7  .  10^^  (297°K) 

-- 

A 

3.3  .  10^^  (297°K) 

-- 

9 

Wentinck 

^2 

1.2  .  10^^ 

-- 

Harteck^^ 

^2 

6.  1  .  10^^ 

-- 

A 

6.  1  •  10^^ 

.  m, 
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Figure  1  The  Shock  Tube  and  Associated  Instrumentation 
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Propagating  into  N 
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into  a  13%  N,-87%  Ne  Mixture 
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Figure  4  Monochromatic  Interferogram  of  a  M  11.  7  Shock  Wave  Propagating 
into  a  15.  5%  Nt-84,  5%  A  Mixture 
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FIGURE  5  INDEX  PROFILE ,  CM 


FIGURE  6  INDEX  PROFILE  ,  CM 
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APPENDIX 


A  computer  program  for  the  IBM  1620  computer  has  been  formulated  to 
calculate  index  of  refraction  profiles  behind  a  normal  shock.  Coding  of  the 
program  is  for  use  with  the  IBM  three-address  interpretive  routine  known  as 
TAISSIT.  The  basic  equations  and  the  logic  used  in  the  program  are  discussed 
in  this  section,  and  a  complete  listing  of  the  program  is  contained  in  a  later 
section,  together  with  detailed  statements  of  the  operations  performed  at  each 
step,  nomenclature  of  intermediate  results,  lists  of  working  storage  addresses 
and  other  information.  The  listing  is  complete  enough  so  that  any  interested 
person  may  use  or  modify  the  program  with  a  minimum  of  effort. 

There  are  three  essentially  separate  functions  performed  in  the  computa¬ 
tion.  First,  the  equilibrium  values  of  temperature  and  degree  of  dissociation 
(and  related  variables)  are  obtained  Second,  the  temperature,  den-sity,  in¬ 
dex  of  refraction  (and  other  quantities)  are  calculated  which  correspond  to  a 
scries  of  predetermined  values  of  dissociation  (o').  Then  an  expression  for 
dx/do'  is  obtained  for  each  value  of  o,  and  integrated  over  r>'  to  obtain  the 
value  of  linear  distance  x  corresponding  to  each  set  of  values  (o',  T,  o,  M,  N). 
Storage  is  provided  so  that  up  to  14  integrations,  each  utilizing  different  sets 
of  temperature  functions  and  efficiency  factors,  can  be  carried  out  during  a 
single  run. 


I.  EQUILIBRIUM  CALCULATION 

A.  In  order  to  start  the  iterative  loop  used  to  find  the  equilibrium  dissoci 
ation  and  temperature,  estimates  of  .ind  T^  .ire  required.  These  are 
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part  of  the  input  data.  The  value  of  oj  is  momentarily  considered  to  be  fixed, 
and  the  closed  T  Subroutine  is  used  to  find  the  corresponding  values  of  T  in 
the  following  way; 


1.  Using  Tq'  andOfe',  evaluate 

-  1  -  a'  Of* 

"  1  +  a'  ^^2  1  +  at 

1  +  y,»  m2 
— - — ^  /?•- 
Vl-Mj 

2.  M*  =  - 

2^'  -  1 

where  yj'  -  Cy^  +  (1-C)y£ 


and  C  =  mole  fraction  of  diatomic  gas  in  mixture. 
(Subscript  i  denotes  inert  gas,  d  diatomic.  ) 


3, 


Vi'm' 
1  +  o' 


1  + 


p»  pjMj/M' 


4.  The  subroutine  then  compares  the  T  calculated  at  step  3  with  the  in¬ 
itial  guess  value  of  T,  T-,  If  the  difference  exceeds  the  tolerance 
stated  in  the  input,  the  calculated  value  of  T  replaces  Tj  and  becomes 
the  guess  value  for  the  next  iteration  starting  again  at  step  1. 
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This  loop  converges  qui'e  rapidly.  The  difference  AT  =  IT  -  Tj|  de¬ 
creases  by  a  factor  of  5  or  more  with  each  iteration,  provided  that  the 
initial  guess  value  of  T  is  within  2000  K  of  the  actual  value. 

When  a  value  of  T  is  found  which  satisfies  the  established  criterion 
0 

(0.2  K  was  used  in  the  work  reported  here),  exit  from  the  subroutine 
occurs  with  the  final  value  of  T,  p  and  M'  available  for  further 
computation. 

B.  The  value  of  T  determined  above  is  used  by  the  Ke  Subroutine  in  a  7- 
point  Stirling  interpolation  formula"^  to  obtain  In  (T)  and  then  (T).  This 
value  is  stored  and  available  for  later  use. 


The  f  Subroutine  then  evaluates  the  expression 
f  =  1  -  — 


p 


(c  -  a')Mw 


mix 


It  can  be  shown  that  this  function  is  a  monotonically  decreasing  func¬ 
tion  of  a',  and  that  f  =  0  is  a  necessary  and  sufficient  condition  for 
equilibrium. 


D.  The  program  tests  the  sign  of  f  and  then  increments  or  decrements 
(as  appropriate)  the  initial  value  of  aj  by  the  amount  Attg  =  0.  01,  calculated 
as  a  new  starting  temperature  from  the  relations 

AT  =  -.5^,  [2D<,/R.t]  . 

Ti  =  To  +  at, 

and  repeats  the  calculation. 


*  Scarborough,  Numerical  Mathematical  Analysis .  3rdEd.  ,  The  Johns  Hopkins 
Press,  Baltimore,  1955.  pp  70-74. 
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E.  When  f  changes  sign,  indicating  that  a* has  passed  the  equilibrium 
value,  Newton's  method  is  applied.  The  values  of  f  immediately  above  and 
below  zero  are  used  to  determine  the  slope  and  the  approximate  value  of  the 
root: 


a'  =  +  6  a' 


6a'=  f, 


a'  -  a'  , 

n _ n.- 1 

n  f  -  f  , 
n  n- 1 


This  6 a* is  then  used  in  the  equation  of  step  D  to  obtain  an  approximate  value 


of  T  corresponding  to  the  root.  The  value  of  ^a^  in  (D)  is  reduced  by  a  factor 
of  10,  and  the  program  returns  to  (A)  and  repeats. 


F.  When  four  significant  digits  of  a'  have  been  obtained  in  this  manner 
the  calculation  is  concluded.  The  final  values  of  a^  snd  T^,  as  well  as  all 
pertinent  intermediate  results,  are  retained  in  memory  for  later  use. 


II.  PROFILE  CALCULATION 


Following  calculation  and  pointing  of  the  equilibrium  values  the  program 
proceeds  to  the  computation  of  the  index  of  refraction  profile. 
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A.  At  the  shock  front  o'  ■  0.  The  starting  value  of  T  is  obtained 
from  the  expression 


>'k  -  '’ic 
’•k  ■ 

This  expression  assumes  complete  translational,  rotational,  and  classical 
vibrational  equilibrium.  At  best  it  gives  only  an  approximation  to  the  value 
of  T  behind  the  front,  but  only  a  crude  approximation  is  needed  to  start  the 
iterative  loop.  The  accuracy  of  the  fina'  answer  does  not  depend  at  all  upon 
the  accuracy  of  the  assumed  starting  value. 

B  The  T  Subroutine  then  evaluates  T  and  p  as  described  in  I.  A. 

C.  Using  the  value  of  calculated  during  the  final  iteration  of  the 
Equilibrium  Calculation  and  the  T  from  (B)  above,  the  f  Subroutine  evaluates 
the  expression  given  in  I  C. 
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D,  The  program  then  calculates  N  and  certain  intermediate  results 
needed  for  the  evaluation  of  intermediate  results  are  listed  in 

the  working  Storage  assignment  list.  The  nomenclature  utilized  is  a  purely 
arbitrary  one;  the  k's  are  quantities  which  remain  constant  throughout  a  run, 
and  the  q's  are  functions  of  T  and  or'as  shown  in  the  listing. 


E,  The  quantity  (^i)  calculated  from  the  relationship: 


dx 

do' 


J_ 

dg* 

dx 


do'  P 

w 

mix 


j^(c.o')^c^kti^  +  2o' (c-a')kd^+  {c-of'Ul-c)kj  J 


where  the 


.2 

4o'  p 


(c-oflW 


TT 


w. 


mix 


k^i's  are  as  indicated  in  Equations  3  of  the  text. 


At  the  time  the  program  was  originally  written,  the  values  Sj  -  2, 

S2  =  Sj  =  1  were  assigned,  the  r  functions  evaluated,  and  the  numerical 
value  combined  with  other  constants.  The  resulting  expression  was 


When  it  became  necessary  to  include  variations  of  Sj  and  S2  =  S3  in  the  com¬ 
putation  this  expression  was  modified  to  the  form 
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q5(D/RT)'"l 


y  = 


‘l30i'  +  (D/RT)  (qiO^’z'  + 
The  effect  of  this  change  is  to  set 


Sj  =  2  +  nj 


82=32=  1  +  Hj  +  n2 


as  far  as  the  exponent  of  (D/RT)  is  concerned,  but  it  does  not  take  into 
count  the  change  in  the  value  of  the  r  function  in  Equation  3.  Thus  the 

A 

merical  value  labelled  C|  in  the  latter  expression  is  the  equivalent  of 


r(3  +  n,) 


r(2)  ' 

with 

similar  relationships  existing  for  the  other  ^'s 

II 

<0 

r(3  +  nj ) 
r(3) 

1 

1 

A 

O2  = 

r(2  +  nj  +  02) 

^2 

r(2) 

A 

0  = 

^'3 

r  (2  +  nj  +  02) 

3 

r(2) 

ac- 


nu- 


This  fact  was  overhmkcd  in  'he  original  computation  and  necessitated 
"ex  post  facto"  modification  of  the  o  values  associated  with  a  given  profile, 
but  otherwise  cavised  no  difficulty.  Anyone  contemplating  use  of  this  program 
in  the  future  should  modify  it  to  remove  this  possibility  of  misunderstanding. 


This  calculation  is  carried  out  in  a  closed  subroutine  ("Special  Program") 
starting  at  Location  200.  In  the  program  listing,  the  program  shown  utilizes 
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1 


a  set  (nj.  n^.  $3)  for  each  profile  desired.  The  v  values  corre- 

spondinptoeachset  are  calculated  and  placed  in  a  storage  file  for  use  in 
forming  the  integrals. 

F.  To  obtain  the  value  of  X  corresponding  to  each  pair  of  a'  and  N  val¬ 
ues,  we  use 

a'n 

f(^)  • 

“n.l 

n 

a„  =  L  Aa„  . 

i  =  o 

This  computation  is  performed  by  the  Integration  Subroutines.  After  the 

first  four  steps  («'  =  . or^  )  carried  out,  Simpson's  Rule  is 

used  to  compute  a^,  and  the  four  point  integration  formula 


Aa 


n 


=  ¥ 


'<Vi 


-  (y 


n-Z 


^n.l) 


IZ 


+  (y 


n- 1 


is  used  to  evaluate  ^^^2’  f*'®” 

®1  =  *2  ■  ^*Z ' 
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initialize  SKiN  TEST 


The  computation  then  proceeds,  using  the  four-point  form  to  compute 
successive  /^a's  until  a  value  of  a*  i®  reached  which  causes  the  function  f  to 
become  negative,  indicating  that  the  equilibrium  value  has  been  exceeded. 
This  causes  the  integration  routine  to  discard  the  invalid  computation,  and 
Simpson's  Rule  is  again  used  to  compute  j  +  where  the  subscript 

f  denotes  the  final  area  strip  in  the  integral.  Since  j  was  obtained  from 
the  last  valid  four -point  integration,  .^a^  can  be  computed.  The  program 
then  returns  to  the  beginning  and  seeks  new  input. 

"BOOT  STRAP  PROGRAM" 


100 

9 

155 

103 

103 

Read  1  word  from  tape. 

101 

9 

580 

103 

103 

Print  the  word. 

102 

7 

661 

000 

103 

HALT  if  SW  1  is  ON. 

Proceed  to  103  when  "START 

is  pushed. 

103 

[o 

000 

000 

ooo] 

Word  read  from  tape. 

104 

7 

700 

000 

100 

Transfer  to  100  and  repeat. 
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SPECIAL  PROGRAM 


200 

6 

063 

221  1 

[ooo] 

Set  exit. 

201 

6 

064 

000 

932 

Set  integration  counter  (882  +  5N) 

202 

6 

063 

217 

886 

Initialize  answer  file  address. 

203 

6 

062 

205 

230 

Initialize  move  order. 

204 

8 

864 

775 

224 

In  Qj  — »  temp.  stor. 

205 

0 

005 

[ooo] 

225 

Move  (nj,  (Py  ^2'  ^3^  wkng 

location. 

206 

5 

225 

224 

000 

-  nj  In  qj 

207 

8 

874 

000 

ooo 

1 

208 

3 

779 

000 

799 

(q^) 

**5  *  n  ^  ^5  • 

(qi) 

209 

3 

226 

224 

000 

n2  X  Inq^. 

210 

8 

874 

000 

223 

"2 

qj  _  temp.  stor. 

211 

3 

784 

228 

222 

qjQ  X  (p2  — »  temp.  stor. 

212 

3 

773 

229 

000 

kj3  X  (pj 

213 

1 

000 

222 

ooo 

+  ‘»10‘P2 

214 

3 

000 

223 

223 

"2 

X  (qj  )  — •  temp,  stor 
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215 

3 

777 

227 

000 

93  X  fPi 

216 

1 

000 

223 

000 

'll  ^  ^*^10  ^2  ^  3^ 

217 

4 

799 

000 

[ooo] 

(q^qj'”^)  [216]  Vj—  an.file 

218 

6 

012 

205 

005 

Increment  move. 

219 

6 

013 

217 

005 

Increment  answer  file  address. 

220 

6 

134 

217 

205 

If  ans  add  £  final  ans  ads,  repeat. 

221 

7 

700 

000 

[ooo] 

Exit. 
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SHOCK  REDACTION  .  PROGRM 


300 

7  661  300  301 

HALT  if  S^V  1  is  ON 

301 

9  155  000  016 

Load  data  from  tape 

302 

6  061  000  000 

Set  IR  1  to  000 

303 

5  062  000  000 

Set  IR  2  to  000 

304 

5  063  000  000 

Set  IR  3  to  000 

305 

7  661  305  306 

HALT  if  SW  1  is  ON 

306 

5  155  307  307 

Head  "Special  Integration"  load  command  from  tape 

307 

[5  155  200  225] 

(Command  read  from  tape) 

300 

7  661  308  309 

HALT  if  SW  1  is  ON 

309 

5  155  311  313 

Read  three  irords  from  tape  (print  contnands) 

310 

5  5O0  800  816 

Print  input  data  ALTERNATE  ENTRY  POINT 

311 

[5  580  211  21  i] 

(Print 

312 

[5  580  221  22  3 

(Print  /^'s) 

313 

[5  580  231  23| 

(Print  ^'s) 

314 

7  661  314  315 

HALT  if  SW  1  is  ON 

315 

3  840  804  750 

*  —  *  '»ikng  storage 

316 

1  000  848  799 

(i-C,)  X  ^  temp  storage 

317 

3  804  849  000 

C,  X  Yg 

310 

T  000  799  751 

(c,  X  n'  )  +•  wkng  storage 

319 

3  804  850  000 

320 

T  000  799  752 

— >  idcng  storage 

321 

T  000  840  799 

1  — >  temp  storage 

322 

3  752  840  000 

4'  -  1 

323 

4  000  799  762 

(K-')  *•  *)  =•  — >  i»kng  storage 

324 

3  841  751  000 

2  - 

325 

4  000  799  761 

(ay')  -r  =.  1^1  trkng  storage 

326 

3  750  847  763 

(j-c,)  X  a. 5“  =  — »•  iiicng  storage 

327 

3  751  801  799 

y, '  X  /*1,  -»•  tenp  storage 

320 

3  000  801  000 

X  M, 

329 

T  000  840  000 

+-  1 

11-59 

irioig  storage 


330 

4 

000 

799 

764 

^  f) 

-r 

— >.  irioig  storage 

331 

3 

801 

801 

760 

/H,  « 

~  — »  Kicng 

stor^e 

332 

2 

751 

840 

000 

y.'- 1 

333 

4 

841 

000 

000 

^  ^  (r.'-i 

) 

334 

T 

000 

760 

765 

mkng  storage 

335 

3 

802 

751 

766 

=  — >  licng 

storage 

336 

3 

800 

801 

767 

=  k-f  — >  Tikng 

storage 

337 

3 

750 

810 

799 

(i-c,)  « m 

urL  — ►  ten5>  storage 

338 

3 

8o4 

805 

000 

C.|  K 

339 

T 

000 

799 

753 

■C<)  /HwL  = 

_ ^  mkng  storage 

340 

4 

851 

000 

798 

•o‘ 

— ►  temp  storage 

341 

3 

000 

799 

000 

X 

(l-C,) 

342 

3 

000 

811 

768 

X 

Kl  =  4,  — > 

mkng  storage 

343 

3 

798 

805 

000 

344 

3 

000 

806 

769 

X  Kj  =5  k.5  ■  * 

mkng  storage 

345 

8 

835 

805 

000 

346 

3 

000 

753 

000 

K 

347 

3 

000 

800 

000 

X 

h 

348 

3 

000 

803 

000 

» 

349 

4 

000 

852 

770 

•  1 
•  1 

(tt  =1  — >  Tikng  storage 

350 

4 

843 

753 

771 

A  ~  n,;j| 

=  4n  — ►  mkng  storage 

351 

3 

841 

809 

759 

Z  X  (U/R) 

— >  mkng  storage 

352 

7 

661 

352 

360 

HALT  if  SW  3 

.  is  ON. 

353 

Q 

000 

000 

000 

354 

5 

000 

000 

000 

355 

0 

000 

000 

000 

356 

0 

000 

000 

000 

NOT  US3D  BY  «10GRAI» 

357 

5 

000 

000 

000 

358 

5 

000 

000 

000 

359 

0 

000 

000 

000 
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SQUIUBRIUIf  CALCULAgON 


360 

5  062  366  370 

Initialize  sign  test 

361 

5  063  366  377 

Initialize  sign  test 

362 

5  002  813  754 

Uove  initial  Tg  to  iricng  loc 

363 

5  000  853  757 

Move  initial  to  AX'^irfcng  loc 

364 

5  069  000  004 

Set  IR  9  to  4  (nuniber  of  sig  digits  in  ot^) 

365 

7  820  367  480 

Call  I  subroutine.  (Returns  ifith  T  in  acc) 

366 

7  760  [000]&0(3 

Sign  test  f  and  transfer 

367 

7  820  368  445 

Call  K.  Subroutine 

8 

368 

7  820  366  545 

Call  f  Subroutine*  return  to  366 

369 

7  600  000  000 

(Not  used) 

370 

5  063  366  381 

HSR3  IF  f  IS  +  CXi  FIRST  PASS 

Re-set  sign  transfer  for  break-out 

371 

7  733  373  372 

HSRS  IF  f  IS  -ON  2nd-l«>h  PASS 

372 

5  580  869  871 

(Optional  print  f*o(*  T  if  SX  3  is  ON) 

373 

5  003  869  865 

Uove  f>  oc*  T  to  storage  location 

374 

0  000  757  000 

Bring  acc 

375 

7  820  365  467 

Call  "Bstinate  T"  subroutine*  return  to  365 

376 

7  600  000  000 

NOT  USED 

377 

5  062  366  381 

HERE  IF  f  IS  —ON  FIRST  PASS 

Re-set  sign  test  for  break-out 

378 

6  063  366  371 

Re-set  sign  test  for  2nd-Nth  pass  negative 

379 

3  855  757  757 

Set  negative  (  x.  -/ 

380 

7  700  000  371 

Go  to  371 

381 

7  733  383  382 

H3Ri  ViIEN  SIGN  OF  f  CHAIDES 

(Optional  print  f*  cc*  T  if  S'ff  3  is  ON) 

382 

5  580  869  871 

383 

6  062  366  370 

Re-inltiallze  sign  test 

384 

6  063  366  377 

Re-ini ti all ze  sign  test 

385 

B  894  757  757 

Set  bd^positive 

386 

5  029  000  001 

Decrement  IR  9 

387 

7  661  387  388 

HALT  if  SW  1  is  ON 

388 

5  509  389  396 

If  IR  9  **0*  Equilibrium  calc  is  finished. 
If  IR  9  ^  0*  decrement  2Wi^and  continue 

389 

4  757  854  757 

10  ->  AOC'' 
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390 

2  869  865  799 

/n-i  — storage 

391 

2  866  870  000 

<  - 

392 

4  000  799  000 

393 

3  000  869  758 

X  =  ^Oi  T*cng  storage 

394 

7  820  365  467 

Call  "Estimate  T"  svdb-routine^  return  to  365 

395 

7  600  000  000 

(Not  used) 

396 

3  872  780  873 

HERE  '.VHBN  FINISHED 

/O'  X  =*  ^  answer  file 

Print  f ,  h  N 

397 

5  58o  869  873 

398 

7  661  398  400 

HALT  if  S-i?  1  is  ON 

399 

7  600  000  000 

(Not  used) 

START  OF  MAIN  PROGRAJi 

400 

6  062  441  555 

Set  integration  Transfer  Sftitch  for  1st  pass 

401 

0  000  857  754 

Initialize  ex.'  to  0 

402 

0  000  816  757 

Initialize 

403 

4  000  841  787 

Aoc*  -f  2.  — ►  itoig  storage 

404 

4  757  842  788 

/>«.*  -i-  5  — ^  wkng  storage 

405 

3  808  808  776 

»  ^ci  --*•  storage 

406 

3  000  856  772 

A  4 -  4,2  — ■"laig  storage 

407 

T  805  810  000 

408 

4  000  810  000 

-f 

409 

8  835  000  799 

^  ^  temp  storage 

410 

3  804  776  797 

X  3>4  — ^  temp  storage 

411 

3  812  812  000 

Di  %  > 

412 

3  000  750  000 

X  (.l-c,) 

413 

T  000  797  000 

-4-  C,-vf- 

414 

3  000  799  000 

415 

3  000  750  000 

X  (/--c.) 

416 

3  000  858  773 

X  2/2”  =  4,3  — w  wkng  stor^e 

417 

7  600  000  000 

No  operation 

418 

0  000  797  785 

Move  a,  into  (c, -‘x’)  location 

419 

3  760  761  799 

A/),*’  *  )(^  - ►  temp  storage 
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420 

2  000  762  798 

tenp  stor£|ge 

421 

3  000  762  000 

422 

T  000  840  000 

1 

423 

4  000  799  000 

-r  4. 

424 

3  000  798  000 

425 

3  000  802  755 

X  T,  7^  — >  Ttaig  storage 

426 

7  820  443  480 

Call  T  Subroutine*  return  to  443 

427 

7  760  428  670 

Sign  test  f.  If  neg*  the  calc  is  invalid*  go  to 
final  integration.  If  pos*  proceed. 

428 

8  874  775  799 

3^9?  (91.)  toJH>  storage 

429 

8  835  871  000 

rr 

430 

3  000  775  000 

431 

4  000  799  000 

-  =  h 

432 

3  000  872  000 

V 

433 

3  000  872  000 

^ 

434 

3  000  774  000 

<  (c,-oc'J 

435 

3  000  869  000 

- 

436 

4  770  000  779 

kfo  ~  i  =^-  — ►  irimg  stori^e 

437 

3  870  772  784 

oc  X  4,2,  f  10  '  >  nloig  storage 

438 

3  785  775  777 

[Cc,  X  storjge 

439 

3  872  780  873 

P»  ^  ^  anseer  file 

440 

7  661  440  441 

HALT  if  SWl  is  ON 

441 

7  820  [64|  200 

Call  "IfultiidLe  Integration  Sub-routine*  then  out 
to  Integration.  The  B  ads  of  this  instruction 

442 

7  600  000  000 

is  the  "Integration  Transfer  Snitch." 

443 

0  000  479  000 

Move  Kq  to  acc 

444 

7  820  427  545 

Call  f  Subroutine*  return  to  427 

K_  INTERPOLATION  SUB-ROUTINE 

445 

6  063  464  000 

Set  sub-routine  exit 

446 

2  000  985  000 

T-T^ 

447 

4  000  986  998 

“  li  =  M.  »•  ten^  storage 

448 

3  000  000  000 

U  X  u.  =  u.^ 

449 

2  000  840  997 

U.^  -  1  tenp  storage 
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450 

5  000  842  996 

(u^-/)  -3  ten?)  storage 

451 

3  998  993  000 

U  * 

452 

T  000  992  000 

453 

3  000  996  996 

A  (u^--4-)  ■=  0^  *■  tenqp  storage 

454 

3  998  991  000 

U  X  Ci, 

455 

T  000  990  000 

+  <i3 

456 

T  000  996  000 

'o-'' 

457 

3  000  997  996 

,2  V  *  L  "  _ 

X  (u  - 1 ;  =  b  — ►  tenqp  storage 

458 

3  998  989  000 

459 

T  000  988  000 

+  C, 

460 

T  000  996  000 

■h  b  —  c 

461 

3  000  998  000 

462 

T  000  987  000 

+  /c  =  .//»  ^  cri 

463 

5  874  000  479 

Bsqp  (In  Ke)  =  Kg  — ^idang  storage 

464 

7  700  000^45] 

£xit  from  siil>-routine 

465 

466 

7  600  000  000 

7  600  000  000 

(Jtot  used) 

■ESnUATS  T"  SU3-S0UTI13 

467 

6  063  475  000 

Set  sub-routine  exit 

468 

Q  000  000  799 

Act*  — ^  tenro  storage 

469 

T  844  870  798 

7  •*-  oil  teiqp  storage 

470 

T  759  871  000 

[aU/R]  +  T; 

471 

3  000  799  000 

A  Aoc' 

472 

4  000  798  000 

■r  Oftoi'i)  =  aT 

473 

2  871  000  755 

^  TTfi  ^  ■otaig  storage 

474 

T  754  799  754 

=*  •»"  *  TJkng  storage 

475 

7  700  000  [423 

iSxit  froD  sub-routine 

476 

477 

478 

479 

0  000  000  000 

0  000  000  000 

S  000  000  000 
[6  000  000  000] 

(Kot  used) 

Working  storage  location.  Kg 

11-64 


T  SUBROUTINE 


480 

5  063  544  000 

Set  sub-routine  exit 

481 

2  8o4  754  774 

C  j  tc*  — >  'vdoig  storage 

482 

T  840  754  799 

1  ^  tenp  storage 

483 

4  766  799  783 

-f  storage 

484 

3  807  754  000 

<»/  ^  oc' 

485 

T  000  804  000 

-H  C, 

486 

J  000  769  000 

^  It, 

487 

T  000  768  780 

+-  4^  =s  — ►  Tilcng  storage 

488 

3  771  754  000 

X 

489 

3  000  754  000 

X  oc' 

490 

4  000  774  778 

T*  (c,  -ot)  =  — ►  nkng  storage 

491 

4  809  755  775 

-t  T  =*  »  Tikng  storage 

492 

T  8o4  754  000 

C,  +  oc^ 

493 

4  804  000  781 

C,  ^  Ce-i  +  x^)  ~  storage 

494 

5  000  755  000 

Uove  T  to  ecc 

495 

9  000  730  720 

Table  Look  Up:  (t) 

496 

2  755  002  799 

T  temp  storage 

497 

2  001  002  000 

T...  -  -‘‘T 

498 

4  799  000  798 

aT  — ^  teiq>  storage 

499 

2  003  004  000 

500 

3  000  798  000 

501 

T  000  004  000 

t-  Aj^CtJ  - 

502 

3  000  774  782 

X  Cc, -«)  =.  ^^>iikng  storage 

503 

0  000  755  000 

lioTe  T  to  acc 

504 

9  000  740  720 

Table  Look  Up:  Cr) 

505 

2  003  004  000 

iSf/Ont-t)  I^H  (T") 

506 

3  000  798  000 

^  '  (V") 

507 

T  000  004  000 

4-  (r« )  =■  (x ) 

508 

3  000  841  000 

/3^  (r)  k  Z 

509 

T  000  775  000 

*  p 
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5tO 

3  000  754  000 

(2/i^cr^  f,)  ^  a' 

511 

T  000  782  000 

h 

512 

3  000  781  000 

^  $r 

513 

T  000  763  799 

+-  ^  storage 

514 

3  000  841  000 

Z  X 

515 

2  000  840  798 

—  1  tenp  storage 

516 

3  000  765  797 

^  Po  ^  storage 

517 

3  764  799  799 

X  ^  temp  storage 

518 

3  000  000  000 

519 

2  000  797  000 

~pO 

520 

8  835  000  000 

521 

2  799  000  000 

^  ‘ 

522 

4  000  798  874 

-r  -  0  ^  XA  — *■  vtsmr  file 

523 

4  767  000  872 

K7  -7-  M  =  ■'**•*“ 

524 

2  764  874  000 

525 

3  000  874  000 

X  M 

526 

3  000  783  871 

X  T  — answer  file 

527 

2  000  755  756 

T  -  ^  storage 

528 

4  809  871  775 

(Jt/R.)  -f  T  *=  |i  — ►  wkng  storage 

529 

7  733  533  532 

Optional  print  T»  aT  if  3  is  OW 

530 

7  600  000  000 

(Not  used) 

531 

7  600  000  000 

532 

9  580  755  756 

Pilot  order  (bypassed  if  SVif  3  is  OFT) 

533 

8  894  756  000 

Set  positire 

534 

2  000  815  000 

I^TI  Terror 

535 

7  661  535  536 

HALT  if  SW  1  is  OW 

536 

7  760  538  541 

Sign  test  and  transfer 

•h  :  Not  finished^  repeat 

537 

7  600  000  000 

•  :  T  calc  finished*  cd.culate  f 

538 

5  000  871  755 

Uore  T  to  Tj^jitial 

'  539 

7  700  000  495 

Repeat 

1  ' 
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she 

7  600  000  000 

(iJot  used) 

541 

0  000  754  870 

Hers  \Then  T  calc  finished. 

Move  oc'to  ans  file 

542 

3  774  776  785 

(c,  -  oc]  A  — »  ukng  storage 

543 

5  000  871  755 

Move  T  to  and  to  acc 

544 

7  661  544  [00^ 

HALT  if  STJ  1  is  ON,  Ihen  exit  from  T  Subroutine 

545 

5  063  549  ^0^ 

f  SUBROUTINE 

Set  exit 

546 

4  872  479  000 

/o  r  i<e. 

547 

5  000  778  000 

548 

5  840  000  869 

^  “  f  answer  file 

549 

7  661  549  [00^ 

HALT  if  S'H  1  is  ON.  Then  exit  from  f  Subroutine 

550 

G  000  000  000 

551 

G  000  000  000 

(Not  used) 

552 

G  000  000  000 

553 

G  000  000  000 

554 

G  000  000  000 

INTEGRATION,  FIRST  PASS 

555 

5  062  441  575 

Set  Transfer  Switch  for  Second  Pass 

556 

5  580  870  874 

Print  oc/ »  T,  ^  ,  N,  Bf 

557 

6  062  561  886 

Initialize  Tj_  address 

558 

5  063  563  889 

Initialize  Yj_  "delqr  line"  address 

559 

G  000  857  790 

Move  0  into  location 

560 

G  000  000  791 

Move  0  into  <1,  location 

561 

G  000  [91  jj  789 

Move  Ypi  from  Hie  to  wkng  location 

562 

5  580  789  791 

/oi  ,  ^ 

563 

G  000  789(915 

Move  Ygj^  to  deli^  line 

564 

6  012  561  005 

Increment  Ygj^  file  address 

565 

5  013  563  005 

Increment  deleor  line  storage  address 

566 

6  124  561  561 

Repeat  if  Y^j^  file  ads  ^  IR  U  (which  contains 

final  address) 

567 

7  661  567  568 

HALT  if  Sti  1  is  ON  • 

568 

0  000  757  000 

Move  acl'  to  acc 

569 

7  820  426  467 

Call  "Estimate  T"  subroutine,  then  return  to 
main  program  at  U26, 
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570 

0  000  000  000 

571 

0  000  000  000 

572 

0  000  000  000 

(Not  used) 

573 

0  000  000  000 

574 

0  000  000  000 

INTEGRATION,  SECOND  PASS 

575 

6  062  441  590 

Set  Transfer  Switch  for  Third  Pass 

576 

0  005  870  880 

Move  T,  N,  H)  to  delegr  location 

577 

6  062  579  886 

Initialize  address 

578 

6  063  579  888 

Initialize  •’delay  line"  address 

579 

0  000  ^li]  &13| 

Move 

580 

6  012  579  005 

Increment  Y^  address 

581 

6  013  579  005 

Increment  "delay  line"  address 

582 

6  124  579  579 

Repeat  if  Y^j^  file  ads  IR  4  (contains  final  ads) 

583 

0  000  757  000 

Move  det'to  acc 

584 

7  661  584  585 

HALT  if  SV.  1  is  ON 

585 

586 

7  820  426  467 

0  000  000  000 

Call  "Estimate  T"  subroutine,  then  return  to  main 
prograo  at  426. 

587 

588 

0  000  000  000 

0  000  000  000 

(Not  used) 

589 

0  000  000  000 

INTEGRATION,  THIRD  PASS 

590 

6  062  441  605 

Set  Transfer  Switch  for  Fourth  Pass  * 

591 

0  005  870  875 

liove  («.^,  T,  ^ ,  N,  M)  to  delay  location 

592 

6  062  594  886 

Initialize  '[2±  eddress 

593 

6  063  594  887 

Initialize  '[2±  "delay  line"  address 

594 

0  000  [91  H  [91  ^ 

k'ove  Y2i 

595 

6  012  594  005 

Increment  Y2i  address 

596 

6  013  594  005 

Increment  Y2i  "del^  line"  address 

597 

6  124  594  594 

Repeat  if  Y2i  ^tle  ads  ^  IR  4  (contains  final  ads) 

598 

7  661  598  599 

HALT  if  S'«  1  is  ON 

599 

0  000  757  000 

Uove  to  acc 

11-68 


600 

7  820 

426 

467 

Call  "Estimate  T"  subroutine,  then  return  to 
main  program  at  h26. 

601 

0  000 

000 

000 

602 

6  000 

000 

000 

(Not  used) 

603 

0  000 

000 

000 

604 

0  000 

000 

000 

INTEGRATION,  FOURTH  PASS 

605 

6  062 

441 

645 

Set  Transfer  S-witch  for  Nth  Pass 

606 

6  062 

609 

886 

Initialize  address  in  del^  line 

607 

5  063 

624  885 

Initialize  new  "delay  line"  address 

608 

9  580 

880 

884 

Print  (ot,  T,  N,  H)  corresponding  to 

609 

0  004  [911]  791 

Move  Y3i>  ...  ,  from  delay  line  to  tikng  loc 

610 

3  843 

793 

000 

4  > 

611 

T  000 

794  000 

612 

T  000 

792 

000 

+  /z 

613 

3  000 

788 

789 

A  (^Vs)  = 

614 

T  794 

791 

799 

)o  +■  >3  *’  ten^j  storage 

615 

T  793 

792 

798 

y,  4-  —*■  tamp  storage 

616 

5  000 

799 

000 

- 

617 

4  000 

845 

000 

12 

618 

T  000 

798 

000 

-»*  (y>  /z) 

619 

3  000 

787 

790 

620 

2  789 

790 

798 

a  2^  —  CL  2  “ 

621 

0  000 

000 

799 

Move  A4.,:  to  print  position 

622 

0  000 

793 

797 

More  Y,i  to  print  position 

623 

9  580 

797 

799 

Print  Y-]  ^ ,  a^  ^ ,  Aaj^^ 

624 

0  005 

789  [91^ 

Uove  (Aaiit Aa2j^), Y^^,  '£2^,  Y^^^  to  delay 

625 

6  012 

609 

005 

Increment  old  delay  line  address 

626 

5  013 

624 

005 

Increment  new  del.8y  line  address 

627 

6  124 

609 

609 

Repeat  if  old  del^  line  ads  £1  IR  4  (final  als) 

628 

7  661 

628 

630 

HALT  if  SW  1  is  ON 

629 

7  600 

000 

000 

(Not  used) 

\ 
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630 

6  062  632  885 

Initialize  del^  line  ads 

631 

9  580  875  879 

Print  (0(,  T,  N>  M)  corresponding  to  Y2 

632 

6  004  (91(3  790 

Move  Y2i  from  delgjr 

633 

5  000  793  789 

Move  Y2i  to  print  position 

634 

9  58o  789  791 

Print  l2i*  ®2i*  ^®2i 

635 

5  012  632  005 

Increment  deify  line  fds 

636 

5  124  632  632 

Repeat  if  delay  line  ads  ^  Hi  4 

637 

0  005  870  875 

Move  (a^  ...  ,  M)  for  Y3 

638 

7  661  638  639 

HALT  if  SW  1  is  ON 

639 

5  000  757  000 

Move  Aft^to  ace 

64o 

7  820  426  467 

Call  "Sstimate  T"  subroutine^  return  to 
main  program  at  426. 

641 

7  600  000  000 

642 

7  600  000  000 

(Not  used) 

643 

7  600  000  000 

644 

7  600  000  000 

INTaOJailON,  N^4  pass 

645 

5  062  648  885 

Initialize  old  delay  line  address 

646 

5  063  659  885 

Initialize  new  delay  line  address 

647 

9  58o  875  879 

Print  (a^,  T»  U,  JI)  corresponding  to 

648 

0  005  [91(3  792 

**0''^®  ®n-l*  •••  »  ^n-3  ^ 

649 

T  796  793  799 

X1-3  /n  — *'  storage 

650 

T  795  794  798 

Xn-i  Xt-i  — >  temp  storage 

651 

2  000  799  000 

(X/t-z  ^  /o-t)  {Xn-3  *  Xn) 

652 

4  000  845  000 

'T  H 

653 

T  000  798  000 

654 

3  000  787  799 

^  (  Xz)  ten^  storage 

655 

T  792  000  791 

^n-z  7*  =  a.„.,  print  position 

656 

5  000  799  792 

Move  L an_ito  print  position 

657 

5  000  794  790 

Move  to  print  position 

658 

9  58o  790  792 

Print 

659 

5  005  791  [91^ 

®n-l,i»  ^®n-l,i*  ^ni*  ^n-l,i*  ^n-2,i  ^ 
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660 

661 

662 

663 

664 

665 

666 

667 

668 

669 

670 

671 

672 

673 

674 

675 

676 

677 

678 

679 

680 
681 
682 
683 


6  012  648  005 

Increment  old  delay  ads 

6  013  659  005 

Increment  new  delay  ads 

6  124  648  648 

Repeat  if  delay  line  ads  IR  4 

6  005  870  875 

Move  (olf  N,  M)  to  delay  storage 

7  661  664  665 

HALT  if  SW  1  is  ON 

6  000  757  000 

tiove  /&<('to  ace 

7  820  426  467 

000  000  000 

Call  "Estimate  T"  subroutine,  return  to 
main  program  at  426. 

5  000  000  000 

(Not  used) 

5  000  000  000 

INTEGRATION.  FINAL  UNE 

6  062  672  885 

Initialize  del^  line  ads 

5  580  875  879 

Print  (  0(.i  T,  ^  ,  N,  M) 

0  005  [93^  790 

Move  a^_p  Y^g  ^  vkng 

3  843  793  000 

^  ^  /»-. 

T  000  794  000 

T  000  792  000 

/« 

3  000  788  000 

2  000  791  791 

(Aa,.,  A«„  )  -  Aa„.,  = 

T  000  790  790 

0  000  792  789 

Move  Yn  to  print  position 

5  580  789  791 

Print 

6  012  672  005 

Increment  delcy  line  als 

5  124  672  672 

Repeat  if  delay  line  ads  £  IH  4 

7  662  683  100 

HALT  if  SW  2  is  ON,  Then  go  to  "Bootstrap" 
for  nesr  data. 
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TABLE 


Location 

T‘‘K 

Location 

720 

5,  000 

730 

721 

6,  000 

731 

111 

7,  000 

732 

723 

H,  000 

733 

724 

9,  000 

734 

725 

10, 000 

735 

726 

11, 000 

736 

111 

12, 000 

737 

^N2 

Location 

% 

4.  249 

740 

2.  557 

4.  306 

741 

2.  622 

4.  353 

742 

2.  703 

4.397 

743 

2,  786 

4.  445 

744 

2.  875 

4.  503 

745 

2.  951 

4.  576 

746 

3.  017 

4.  666 

747 

3.  072 
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constants  used  by  the  program 


Location 


Quantitv 


.<■'4  0 

1.  0 

841 

2.  0 

842 

3.  0 

843 

4  0 

844 

9.  0 

84  5 

12,  0 

846 

1.  50 

847 

2.  50 

848 

y.  =  1 . 6666667 

849 

1.40 

850 

1.2857143 

851 

10^ 

852 

r  1 

No  RyJ 

853 

=  0  01 

854 

10  0 

855 

-1 

856 

^yfV  =  9. 797960 

857 

0 

858 

2  VT=  2.  82842- 

9. 73393  X  10 


27 
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WORKING  STORAGE  ASSIGNMENTS 


Location 

Quantity 

750 

(1  -  c) 

751 

yj  =  c  +  (1  -  c)  y. 

752 

=  c  Vi^  T  (1  -  c)  y. 

753 

M 

w 

mix 

754 

a  (working  value) 

755 

■^initial'  ('^o’^king  value) 

756 

*  ^initial) 

757 

Att'  (working  value) 

7  58 

f)  a'  (from  Equilibrium  Calculation) 

759 

2(D^/R) 

760 

Mi^ 

761 

^1  = 

762 

k2  =  (rk'  •  /  ^’'k 

763 

k3  =  2.  5  (1  -  c) 

764 

^4  =  (1  +  yj  M  j^)  yj  Mj 
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Location 


Quantity 


765 

766 

767 

768 

769 

770 

771 

772 

773 

774 

775 

776 

777 

778 

779 


+  2/  (y;  -  1) 

^^6  =  W 

^7  =  pjMj 

kjj  =  10^  Mw.  (1  .  c)  K.  /Mw  . 

w  1  1  mix 

>'U-  K 

k, ,  =  4  /  Mw  . 

11  mix 

k,2  -  4  d/ 

k,3  .  (1  -  c)  2  Dave^  ['“w^ 

(C  -  a' ) 

qj  (D^/R)/T 

93  =  (c  -  a')  qj 

2 

94  =  ''1 1  Of'  /  (c  -  O') 

'’S  =  ^10  /  ['^2  (c-a')  f] 
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1/2 


I 


Location 

Ouantitv 

780 

^6  =  ^8  N 

781 

=  c  /  (c  +  a') 

782 

qg  =  (C--&M 

783 

qq  =  kg/(l  +  O'*) 

784 

^110=  ^12 

785 

(c  -  Of' )  d/ 

786 

— 

787 

Ao’  /2. 

788 

/N«’/3 
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AUTHOR  SUMCT  OASSMCAnON  MO. 

I.  S.  Scala  Shock  Wave*  R62SD032 

II.  B.  Cary  Non-equilibrium  Ca»  March  30.  196; 

rmf 

HIGH  ALTITUDE  SHOCK  WAVE  STRUCTURE, 
_ Part*  I  and  U  _ _ _ _ 

AatTRACT 

I.  A  theoretical  model  is  developed  for  the  Alrixrlure  of 
a  shock  wave  in  a  diatomic  i;«»  includini?  rotational  and 
vibrational  relaxation  phenomena. 

II.  Experiments  are  carried  out  In  a  COfttbuRtion  driven 
shock  tube.  Experimental  data  for  %*ibranonal  relaxation 
times  and  dissociation  rate*  are  obta^ed  |or 

aov.oAH  I  I 


None 


i«ien*e  Elrctroric* 
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